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Abstract 

GSH is a continuum mechanical theory constructed to qualitatively account for a broad range of 
granular phenomena. To probe and demonstrate its width, simple solutions of GSH are related to 
granular phenomena and constitutive models, including (i) for vanishing shear rates: static stress 
distribution and propagation of elastic waves; (ii) at slow rates: critical state, shear band, the 
models of hypoplasticity and barodesy; (iii) at higher rates: the MIDI-model, rapid dense flow in 
the Bagnold regime. A unified, densely correlated understanding of granular physics emerges as a 
result of these phenomena ordered and explained employing a single framework. 

PACS numbers: 45.70.n, 81.40.Lni, 83.60.La, 46.05. +b 
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I. INTRODUCTION 

A. Engineering Models 

Being a subject of practical importance, elasto-plastic deformation of dense granu 
media has been under the focus of engineering research for many decades if not centuries 
6|]. The state of engineering theories, however, is confusing for physicists: Innumerable 
continuum mechanical models compete, employing strikingly different expressions. In a 
recent book on soil mechanics, phrases such as morass of equations and jungle of data were 
employed as metaphors [sl- Moreover, this competition is among theories applicable only 
to the slow shear rates of elasto-plastic deformation, while rapid dense flow (such as heap 
flow, mud slide or avalanches) is taken to obey yet rather different equations As a result, 
most engineers believe it is illusory to look for a unifled theory capable of simultaneously 
accounting for granular phenomena of arbitrary shear rates: static stress distribution at 
vanishing rates, the rate-independent, elasto-plastic motion at slow rates, and the rapid 
dense flow at high rates. 

On the other hand, although many theories achieve considerable realism when conflned 
to the effects they were constructed for, they are in essence clever renditions of complex 
data, not a reflection of the underlying physics. We therefore thought it may be worthwhile 
to try out a new starting point, by focusing on the physics, while leaving the rich and subtle 
granular phenomenology aside while constructing the theory. The hope was to arrive at a 
theory that, though not necessarily detailed or realistic in every aspect, is widely applicable 
over the compete range of shear rates, flrmly based in physics, and affords a well founded, 
transparent understanding. A tried and true method to achieve this aim is to construct the 
so-called hydrodynamic theory that physicists take to mean the long-wave-length, continuum 
theory of condensed systems - in contrast to engineers, who use the word as a synonym for 
the Navier-Stokes equations. 



B. GSH — A Hydrodynamic Theory 

The hydrodynamic formahsm was pioneered by Landau js^ 



and Khalatnikov jol in the 

context of superfluid hehum, and introduced to complex fluids by de Gennes [l^. Its two 



crucial points are: First, the input in physics that serves mainly to specify the complete 
set of state variables; and second, the simultaneous consideration of energy and momentum 
conservation, complete with their respective fluxes. Combining this with thermodynamic 
considerations, one finds many more constraints, and far less liberty, than the engineering 
approach to constitutive relations, when constructing expressions for the stress. This is 
an advantage especially for granular media: If a choice of random phenomena is rendered 
correctly by a hydrodynamic theory, chances are that the rest is also adequately accounted 
for - because the theory complies with all the mentioned general principles and is founded 
primarily on insights into the basic physics of a system, not a subset of experimental data. 
Hydrodynamic theories [12I have been derived for many condensed systems, including 



liquid crystals 



magnetism 29|- 
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19|, superfluid ^He |20l-l25l|. superconductors |26l-l28|. macroscopic electro 



32| . ferrofluids |33l-l41|. and polymers j42 



45| . We contend that constructing 



a hydrodynamic theory is both useful and possible for granular media: Useful, because it 
should help to illuminate and order their complex behavior; possible, because a significant 
portion has already been done. We call it GSH, for "granular hydrodynamic theory." It 
divides granular behavior into three regimes, with the jiggling of the grains - quantified as 
the granular temperature Tg - serving as a switch: 

• At vanishing shear rates, grains hardly jiggle, and Tg 0. Static stress distribution 
and the propagation of elastic waves are phenomena of this regime. We call it quasi- 
elastic because the stress stems from deformed grains and is elastic in origin. 

• At slow rates, the jiggling increases and Tg is slightly elevated. Although the stress is 
still predominantly elastic, it may now relax: When the grains loose contact with one 
another briefly, both granular deformation and the associated stress will decrease. This 
is the hypoplastic regime where the hypoplastic model jj] and other rate-independent 
constitutive relations are valid. Typical phenomena are the critical state [ij and in- 
cremental nonlinearity, or the strikingly different loading and unloading curves. 



At high shear rates, we have the rapid dense flow behavior covered by the MIDI 
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model 46| and Bagnold flow. The jiggling is so strong that it exerts a pressure, and 
viscosities are important. They compete with the elastic stress, becoming dominating 
at very high rates. 



C. Validity of General Principles 

There are many arguments in the literature contending that granular media, being unique, 
violate general principles, hence do not have a hydrodynamic theory (as flrst conjectured 



by Kadanoff 



47j). At closer scrutiny, none of these arguments is watertight. Four of which 



(rendered in italic) are as briefly refuted here as is appropriate for an introduction. 

1. "The energy is not conserved in granular media." Although the kinetic energy of the 
grains is not conserved, the total energy is, which includes the heat in the grains. 

2. "Fluctuation- dissipation theorem is not valid in granular media. " There are two ver- 
sions of it, one in terms of the granular temperature Tg, the other in terms of the 
true temperature T. The latter is a general principle. It always holds and is equally 
apphcable to a block of copper and a pile of sand, quantifying how much, eg., the 
volume of each fluctuates. The former is an imperfect analogy, not a general principle. 

3. "The Onsager Relation does not hold in sand, because the underlying microscopic 
dynamics, inelastic scattering, is irreversible. " The true microscopic dynamics in sand 
is, as everywhere else, the reversible Schrodinger equation for the constituent atoms. 

4. A sand pile has much more gravitational energy than a monolayer of grains. Only the 
latter is in equilibrium, the minimal energy state. The former, being "jammed" and 
prevented to reach the former, is too far off equilibrium for thermodynamics to hold. 
Similar to two chambers of air separated by a stuck piston, a pile of sand at rest is in 
fact in equilibrium. The air is in equilibrium because all its many degrees of freedom 
are except one: the position of the piston that upholds a constraint on the volume 
of the two subsystems. In a macroscopic body, all elastic degrees of freedom are in 
equilibrium if the force balance holds, implying the sum of gravitational and elastic 
energy is minimal. Two elastic bodies, one on top of another, are also in equilibrium 



if the sum of their energy is minimal - though there is the constraint that the upper 
body must not shde with respect to the lower one. A sand pile is many little elastic 
bodies on top of one another. If they are constrained to stay put, and their total 
energy is minimal, the pile is in equilibrium. 

D. Two- Stage Irreversibility 

To derive the hydrodynamic theory for granular media, one needs the input of what the 
essence of granular physics is. Our working hypothesis is that it is encapsulated by two 
notions: two-stage irreversibility and variable transient elasticity. The first is related to 
the three spatial scales of any granular media: (a) the macroscopic, (b) the mesoscopic, 
intergranular, and (c) the microscopic, inner granular. Dividing all degrees of freedom into 
these three categories, we treat those of (a) differently from (b,c). Macroscopic degrees 
of freedom, such as the slowly varying stress or flow fields, are specified and employed as 
explicit state variables, but intergranular and inner granular degrees are treated summarily: 
Instead of being specified, only their contribution to the energy is considered and taken, 
respectively, as granular and true heat. So we do not account for the motion of a jiggling 
grain, only include its strongly fluctuating kinetic and elastic energy as contributions to the 
granular heat, characterized by the granular entropy Sg and temperature Tg. Analogously, 
a phonon, or any elastic vibration within the grain, are taken as part of true heat, part 
of s and T. There are only a handful of macroscopic degrees of freedom (a), innumerable 
intergranular ones (b), and yet many orders of magnitude more inner granular ones (c). So 
the statistical tendency to equally distribute the energy among all degrees of freedom implies 
that the energy decays from (a) to (b,c), and from (b) to (c), but never (or hardly ever) 
backwards. This is what we call two-stage irreversibility, see Fig[T] 

The system is in equilibrium only if the true entropy is maximal. Maximal granular 
entropy (given in the special cases of the Edward entropy by counting the number of ways 
to stably build a sand pile, see Sec I VI Cll would characterize equilibrium, only if there were no 
energy decay from (a,b) to (c), or when it is slow enough to be neglected. As the ubiquitous 
inelasticity of granular collisions clearly demonstrate, this is never the case. 

A division into three scales works well when they are clearly separated, when the system 
is much larger than the grains. This is indeed a problem with granular media, though 

7 



dissipative processes 
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FIG. 1. Two-stage irreversibility. Dissipative processeses produce either granular entropy Sg, or 
directly thermal entropy S. Eventually, Sg is also converted to S. 

one of accuracy, not viability. Scale separation is usually better satisfied in engineering 
experiments than in some of physicists. Using glass or steel beads, typically larger than 
sand grains, aggravates the problem. Same is true of 2D experiments employing less and 
larger disks. On the other hand, when there is too little room for spacial averaging, one can 
average over time and runs that also get rid of fluctuations not contained in a continuum 
theory. Moreover, one can go to higher order gradient terms, as we shall do in Sec IVII Ci 
to capture qualitatively what happens at small length scales - eg. in shear band or when 
clogging occurs. 



E. Variable Transient Elasticity 

Our second notion, variable transient elasticity, addresses granular plasticity. The free 
surface of a granular system at rest is frequently tilted. When perturbed, when the grains 
jiggle and Tg ^ 0, the tilted surface will decay and become horizontal. The stronger the 
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grains jiggle and slide, the faster the decay is. We take this as indicative of a system that 
is elastic for Tg = 0, turning transiently elastic for Tg ^ 0, with a stress relaxation rate that 
grows with Tg. A relaxing stress is typical of any viscous-elastic system such as polymers. 
The unique circumstance here is that the relaxation rate is not a material constant, but a 
function of the state variable Tg. As we shall see, it is this dynamically controlled, variable 
transient elasticity - a simple fact at heart - that underlies the complex behavior of granular 
plasticity. Realizing it yields a most economic way to capture granular rheology. 

Employing a strain field rather than the stress as a state variable usually yields a simpler 
description, because the former is in essence a geometric quantity, while the latter contains 
material parameters such as the stiffness constant. Yet one cannot use the standard strain 
field eij as a granular state variable, because the relation between stress and ejj lacks unique- 
ness when the system is plastic. Engineering theories frequently divide the strain into two 
fields, elastic and plastic, eij = Uij + with the first accounting for the reversible and 
second for the irreversible part. They then employ e^j and e[f as two independent strain 



fields to account for granular plasticity 



48 



49|. 



We believe that, on the contrary, the elastic strain Uij is the sole state variable, as there is 
a unique relation between the elastic stress Uij and Uij, if the latter is appropriately defined 
via the elastic energy: Shearing a granular system, a portion of the strain goes into deforming 
the grains individually, changing their elastic energy. The rest of the strain is spent sliding 
and rolling the grains. Taking Uij as the portion that changes the energy and deforms the 
grains, the energy w is by definition a function of Uij alone. And since an elastic stress Hij 
only exists when the grains are deformed, it is also a function of Uij. Therefore, we employ 
Uij as the sole state variable, and discard both eij and el^K Doing so preserves many useful 
features of elasticity, especially the (so-called hyper-elastic) relation. 



-dw{uij)/duij. (1) 



This is derived in 



50| but easy to understand via an analogy. Driving up a snowy hill slowly, 
the car wheels will grip the ground part of the time, slipping otherwise. (We assume a slowly 
turning wheel and quickly changing, intermittent stick-slip behavior.) When the wheels do 
grip, the car moves upward and its gravitational energy ly^*"™ is increased. If we divide the 
wheel's rotation into a gripping (e) and a slipping (p) portion, = 9^^^ + 6^^\ we know we 
may ignore 9^p\ and compute the torque on the wheel as dw^'^"''" / d6^^\ if the wheel turns 
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sufficiently slowly, same as Eq ([T]). How much the wheel turns or slips, how large 6 or 6'^^^ 
are, is irrelevant for the torque. 

The only way to find out whether our two hypotheses are appropriate and complete, is to 
derived the theory and compare its ramifications with experiments. The theory has already 



been derived, see 
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5l| 5 and is called GSH. It is briefiy repeated and presented in the next 



chapter. The second step, finding the ramifications of GSH is a more lengthy process, in the 
midst of which we are. And this manuscript is an overview on the work done and planned. 
As mentioned, GSH has three rate regimes, given by: 



• The quasi-elastic regime of vanishing shear rates, possibly below 10~^ s~^, with a 
quadratically vanishing granular temperature, Tg ~ 7^ — > 0. The stress (Tij = TTij is 
purely elastic, as given by Eq ([1]). 

• The hypoplastic regime of low shear rates, possibly between 10~^ and 1/s, where the 
engineering theory of hypoplastic model [4^ holds. The stress (7^ = (1 — a)7iij is still 
elastic, but softer by a factor typically between 0.2 and 0.3. Granular temperature is 
more elevated, allowing stress relaxation. With Tg ~ 7, we have rate-independence. 
This regime is frequently hailed as the quasi-static one, because it seems slow, is 
rate-independent, and because the even slower quasi-elastic regime is (as we shall 
see) hard to observe. We note that the hypoplastic regime, being characterized by 
stress relaxation, is dissipative. It therefore cannot possibly be quasi-static, implying 
a consecutive visit of neighboring equilibrium states. 

• The rapid flow regime, for shear rates well above 1 s~^. We still have Tg ~ 7, but it is 
no longer small. Therefore, the T^-generated, seismic pressure Pt ^ Tg 7^ and the 
viscous shear stress as ~ Tg'j ~ 7^ become significant and compete with the elastic 
contribution vTjj. This is where the MIDI model and Bagnold fiow hold. As both the 
pressure and the shear stress may be written as ci + 627^, where ei is the elastic, and 
62 the seismic or viscous, contributions, we have the purely quadratic dependence of 
the Bagnold fiow for 627^ ^ ei, and hypoplastic rate- independence for 627^ ^ Ci. 
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II. BRIEF PRESENTATION OF GSH 



A. Complete Set of State Variables 

In accordance to our understanding of granular media's basic physics, its state variables 
are: the granular entropy Sg and the elastic strain Uij, in addition to the usual variables: 
the density p, the momentum density pvi, the true entropy s. Denoting the energy density 
(in the rest frame, Vi = 0) as w = w{p, s, Sg, Uij), we define the conjugate variables as: 

^=~d^' ^=97' '^d^g' ""''^'d^j' ^ ^ 

where p is the chemical potential, T the temperature, Tg the granular temperature, and TTj^ 
the elastic stress. These are given once the energy w is. Next, in Sec III B 1\ equilibrium 
conditions will be derived, formerly, in terms of the energy and its conjugate variables, 
whatever w is. Then, in Sec III B 2\ an example for w will be given, and the conjugate 
variables calculated - with the help of which the equilibrium conditions are rendered explicit. 

A complete set of state variables is one that determines a unique macroscopic state of 
the system. If a set is given, there is no room for ambiguity, for "history-" or "preparation- 
dependence." Conversely, any such dependence indicates that the set is incomplete, see eg. 
the discussion in Sec IVIBI In the hydrodynamic approach, a physical quantity is a state 
variable if (and only if) the energy w depends on it. We assume the above set is complete. 

Having specified the thermodynamic energy w as a function of relaxing variables such as 
Sg, we employ in effect a generalized notion of equilibrium, and treat a state with a finite Sg as 
being in quasi- equilibrium. From a statistical mechanical point of view, this is a constrained 
equilibrium, because we are considering only those micro-states that are compatible with 
the given value of Sg. (An example for such a relaxing thermodynamic variable is the 
magnitude of an order parameter in a Ginzburg-Landau theory, say the superfiuid density 
Ps, cf. [9].) Such a variable needs to be macroscopically slow, so microscopic variables have 
time to adjust to its value. Since Sg typically varies on the scale of 0.1-1 ms in the dense 
limit, much slower than any microscopic time scales, Sg is a valid macroscopic and quasi- 
thermodynamic variable. [The notion of quasi-equilibrium also holds as local equilibrium, 
implying w[r, t] = w[p{r,t), s{r,t), Sg{r,t),Uij{r,t)]. Analogous equations hold for all the 
conjugate variables, see Eqs ([2]).] 
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B. Statics 



1. Equilibrium Conditions 

Requiring maximal entropy J s d^r with appropriate constraints (such as given energy 
J wd^r and mass J pd^r), one obtains the equihbrium conditions for the state variables 
in terms of their conjugate variables. In granular media, remarkably, this universally valid 



procedure leads to two distinct sets of equilibrium conditions, the solid and the fluid one 
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52| . Maximizing the entropy, we first obtain the condition of uniform true temperature 
VjT = 0, and the requirement that the granular temperature vanishes, Tg = 0. Usually, 
Tg vanishes quickly, and if it does, the density is not independent from the elastic strain, 
dp/p = —dufj. They share a common condition that we identify as the solid one, 

Vi{7iij + PrSij) = pgi, (3) 
Pt = -d{wV)/dV 0, (4) 

where gj is the gravitational constant, Hij the elastic stress, Pt the usual expression for the 
fluid pressure, and V the volume. (The derivative is taken at constant M = pV.) We may 
equivalently calculate Pt as —d{w/p)/d{l/p) = p^d{w / p) /dp, holding constant s/p,Sg/p. 
With the energy expression w of the next Sec III B 21 Py^ ~ Tg is the pressure exerted by 



jiggling grains. We therefore call it the seismic pressure [52|- Clearly, equilibrium condition 
Eq ([3]), expressing force balance, is logically the result of maximal true entropy. 

If Tg is kept finite by external perturbations, the system may further increase its entropy 
by independently varying p and Wjj, to arrive at the fluid equilibrium. It is characterized by 
two conditions, the first with respect to Uij, and the second with respect to p: 



TT. 



0, V,Pt = P&. (5) 



The first condition requires shear stresses to vanish in equilibrium, and free surfaces to be 
horizontal. The second condition is that governing reversible compaction, a phenomenon 
one arrives at after keeping Tg finite (such as by tapping) for a long time, see Sec IVll 

2. Granular Energy 

Interested in hard grains that are slightly excited, implying small Uij,Sg, we look for 
the respective lowest order terms in the energy. (As we are not, at present, interested in 
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thermal effects such as thermal expansion, the energy's dependence on the true entropy is 
not discussed.) Denoting A = — Pa = 7r«/3, = u*jU*j, txI = T^tjT^ij, where u*j,7!'*j are 
the respective traceless tensors, we take the energy to be 

w = wt + (6) 

WT = sy{2ph), T, = s,/ipb), (7) 

WA = VA(2i3AV5 + Aul), (8) 

TTij = VA{BA + Aul/2A)6ij - 2AVAu*j, (9) 

Pa = VA{BA + Aul/2A), vr, = -2AVAus. (10) 



Note Uij and tt^j are colinear and have the same principal axes. The contribution wt is an 
expansion in Sg. The quadratic term is the lowest order one because we require Sg,Tg = 
to be a minimum of the energy. (This is the same argument as in a Ginzburg-Landau 
expansion, though without the fourth order term or the phase transition.) As it will turn 
out, see Eq f|T2|) below, this lowest order term is in fact sufficient to account for fast dense 
flow and the gaseous state. In this sense, Sg and Tg are always small. 

Next, we compare Tg to the gaseous granular temperature Tg, defined as 2/3 of the kinetic 
energy per particle, see jssj]. Being a general expression for granular heat, wt includes the 
quickly fluctuating part of both the kinetic and elastic energy. But in the dilute limit, when 
the elastic contribution may be neglected, one can take wt as equal to 3pTG/2{m) (with 
(m) the average mass of a grain), and identify 

3TG/{m) = bTl (11) 

Fixing the density-dependence of the coefficient b immediately yields an expression for 
the seismic pressure Pt = —d{w/p)/d{l/p). [There is also a contribution from wa ~ A^'^, 
because A, B depend on the density, see Eq (fT6|) . It is neglected because it is always much 
smaller than Pa ~ A^'^ for small A.] We take 

pY p^abT^ , ^ 

Pcv) 2(pcp-p) 
with both 6o and a being positive numbers. Given Eq ( ITT]) (noting the density dependence 



of 6), this is essentially the famihar pressure expression ~ Tg/{Pcp — p), see eg. [54 1- 

The second term wa of Eq ([7]), with A,B>0, is the elastic contribution. Its order of 
2.5 is important for many granular features, especially stress-induced anisotropy (see below) 
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and the convexity transition, discussed in Sec III C[ The associated stress expression vTjj has 
been vahdated for the following circumstances, achieving good to satisfactory agreement: 



Static stress distribution in three classic geometries: silo, sand pile, point load on a 



granular sheet, calculated using the equilibrium condition, Eq ([2]), see |56l. |57|. 



Small-amplitude stress-strain relation, see j58|, |59| 



Anisotropic propagation of elastic waves, see [60|, IM] 



An explanation of "stress-induced anisotropy": In linear elasticity w ~ , we have constant 
second derivatives d'^w/duj., and the velocity of a elastic wave ~ a/ d'^w/dul does not depend 
on the strain, or equivalently, the stress. For any exponent other than 2, the velocity depends 
on the stress, and is anisotropic if the stress is. 

Note that the energy w = wt + vanishes when the grains are neither deformed nor 
jiggling: w — > for Sg, Uij — )■ 0. This implies a lack of interaction among the grains. If there 
were any, there would be a third term in w that is a function of p alone. 



C. Yield Surfaces 



In a space spanned by stress components and the density, there is a surface that divides 
two regions in any granular media, one in which the grains necessarily move, another in 
which they may be at rest. This surface is usually referred to as the yield surface. Aiming 
to make its definition more precise, we take the yield surface to be the divide between two 
regions, one in which elastic solutions are stable, and another in which they are not. Clearly, 
the medium may be at rest for a given stress only if an appropriate elastic solution is stable. 
Since the elastic energy of any solution satisfying the equilibrium condition Eq (jS]) is an 
extremum, the energy is convex and minimal in the stable region, concave and maximal in 
the unstable one — in which infinitesimal perturbations suffice to destroy the solution. 



1. The Coulomb Yield Surface 



The elastic energy of Eq ([8]) is convex only for 

Us/ A < sJlBjA or tt./Pa < \/lAIB, 
14 



(13) 



turning concave if the condition is violated. The second constraint may be derived by 
rewriting Eq f lTOj) as 

which shows Pa/t^s = ^/B/2A is minimal for Ms/ A = y^2B/A. This corroborates the 
behavior that no granular system stays static if the shear stress is too large for given pressure. 
We take B/A to be density independent and approximately 

B/A ^ 5/3. (15) 

Therefore, we only need to specify the density dependence of B{p), taking it as 

B = Bo[ip-p)/ip.,-p)r^ (16) 
p = (20p,p - npcp)/9, 

with ,Bo > a material constant. This expression accomplishes three things at once: 

• The energy is concave for any density smaller than the random loose one pep, implying 
no elastic solution exists there. 

• The energy is convex between the random loose density pgp and the random close 
one Pep, ensuring the stability of any elastic solutions in this region. In addition, the 



density dependence of sound velocities as measured by Harding and Richart 62| is 
well rendered by \/B. 

The elastic energy diverges, slowly, at pcp, approximating the observation that the 
system becomes orders of magnitude stiffer there. 



2. Yield Stress versus the Critical State 



A widespread confusion is addressed in this section. The yield surface of Eq f|T3|) defines a 
yield shear stress for a given pressure. Many textbooks identify this stress with the highest 
stress achieved in an approach to the critical state, and draw conclusions based on this 
identification. Their justification is that the approach is typically executed at low enough 
shear rates to be considered quasi-static. We contend that a true quasi-static motion is one 
that visits a series of static, equihbrium states, with — )■ 0. This happens, as mentioned 
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above, only during quasi- elastic motion, see also Sec lIII Al The rate- independent, hypoplastic 
motion, taking place during an approach to the critical state, produces an elevated Tg and is 
strongly dissipative. Therefore, the instability discussed here and the critical state discussed 
in Sec IIVBI are two distinct concepts, static versus dynamic. 

The first is a convexity transition of the elastic energy, to be probed by quasi-elastic 
motion at vanishing Tg. The second is a stationary solution of the evolution equation for 
the elastic strain Uij, and is comparable to the stationary solution of any diffusion equation. 
The two shear stresses are frequently similar in magnitude, but the yield stress given by 
Eq ( |T3l) needs to be larger than the highest shear stress achieved during the approach to 
the critical state, see Fig [5] below. Otherwise, the system will abandon the approach and 
develop shear bands instead, considered in Sec IVIIi 




FIG. 2. Granular yield surfaces for Tg = 0, as a function of the pressure P, shear stress ag = 
TTg, and the void ratio e, as calculated from the convexity transition of the energy: w/Bq = 
^/A (|A2 + lul) - (Pi A3 + Psn^A + V^ui) , with ? ^ | = 5/3, V, = 1, ^2,2^3 = 2, 
and = 0.67 Pep (implying pip = 0.85pcp). The (a,b) are at ag = 0, where the inset has a 
logarithmic scale, while the curves of (c) are at the indicated densities. The dashed straight lines 
in (b,c) are, respectively, the engineer formula e = — klnP and the Coulomb yield line. The 
curves of (d) are the same as in (a,c), though now in the space spanned by e,P, 0"^,. So this gives 
the boundary surface of static states. 
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3. The Virgin Consolidation Surface 



As depicted in Fig. |21 granular media possess more yield surfaces in the space spanned 
by the pressure Pa, shear stress as = vr^, and the void ratio e, where e = 1/0 — 1 {(p = p/pg 
is the packing fraction, and pg the bulk density of the grains). First, for given e, there is 
a maximal pressure that a granular system can sustain before it collapses, implying a yield 
surface as depicted in (a) of Fig[2l This is a boundary that sand at rest will not cross when 
compressed. Instead, it will collapse, becoming more compact, with a smaller e, coming 
to rest at a point close to the curve, never above it. In soild mechanics textbooks, the 
boundary is frequently referred to as the virgin consolidation line. (The inset, (b) of Fig|21 
has a logarithmic scale. It serves to demonstrate that the standard formula e = — klnP 
from engineering textbooks do not go to pip and pcp, for P — > and oo, respectively.) 

This qualitative behavior is easily accounted for by higher-order strain terms. Taking 
Vi,V2,T>3 > as functions of the density, we choose 

- Bo{V,A^ + V.Aul + V,ut), (17) 

to be added to w^, Eq ([8]). Consider first = 0. If A is large enough, the term —ViA^, 
with a negative second derivative, will work against BA"^'^ and turn w/^ concave. The value 
of A at which this happens, call it A^ is given by y/A^ = 5B{p)/8Vi{p). As B diverges at 
Pep, so does Ac. If Ac(p) = for p = p^, T>i{p) will have to diverge there. 

Next consider ^ 0. If V2,'D3 = 0, the yield lines in the space spanned by Pa,tTs for 
given density would be vertical lines. The presence of —V2AUI and —V^uj reduce the value 
of A (or Pa) for growing Ug (or tt^), bending the lines to the left. Although qualitative 
figures of these curves - frequently referred to as "caps" - abound in textbooks , we 
did not find enough granular data, especially not a mathematical expression that we could 
have used to fix Vi,V2,T>^. 



4. More Yield Surfaces 



The yield surface of Eq (1131) is usually referred to as the Drucker-Prager approximation 
of the Coulomb yield surface. The actual Coulomb law is different. And there are more 
yield laws, such as Lade-Duncan [63,] or Matsuoka-Nakai |6J]. Engineers choose among 
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them depending on the system, personal preferences and experiences, apparently without a 



commonly accepted rule. By including the third strain invariant = u*m*^u*^^ into Eq (jH]) 



with A,B,C > 0, it is in fact possible to account for many of them at the same time. 
(Note the new term is of the lowest order, 2.5.) Depending on how large C/B is, the 
convexity transition takes place at yield surfaces that are numerically indistinguishable from 
the respective yield law. Because a single expression is employed, and because intermediate 
yield laws are also possible, this is a simplifying and unifying step. Evidence that these 
laws may be unified using Eq ( fTSl) . and that Uij, iiij retain their colinearity, is given in [65|. 
Defining the friction angle as = arcsin A/3/(6P^/7r^ — 1), the Coulomb, Drucker-Prager, 
and Lade-Duncan yield laws are given respectively as 

TTg-VTi . TVs y/6 Sm 

, = smip, — = ^= , (19) 

+ TTi Pa VS + sinV 

27(1 — sin (p) cos^ if 

Pi ^ (3-sin(^)3 ■ 

D. Dynamics 

1 . Structure of the Dynamics 

Next, we specify the evolution equations for the state variables. The equation for the 



elastic strain, assuming both uniform, is |50|] 



dtUij - Vij + aijkeVM = -{.\juTg) uu (20) 

(where Vij = \iSiVj + ^ jVi) is the shear rate, v*^ its traceless part, and vl = v*jV*j). If Tg is 
finite, grains jiggle and briefly lose contact with one another, during which their deformation 
is partially lost. Macroscopically, this shows up as a relaxation of Uij, with a rate that grows 
with Tg, and vanishes for Tg = 0. So the lowest order term in an T^-expansion is XijkeTg. 
With the elastic energy a convex function, the (negative) elastic stress —i^ij = dw/duij is 
a monotonically increasing function of Uij. Therefore, —iTij^Uij decrease at the same time. 
And Eq ( 120|) accounts for the stress relaxation discussed in the introduction. 
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The Onsager coefficient aijke is an off-diagonal element. Dividing Uij into A = —uu, u*j, 
and specifying the matrices a^ke, \jM with two elements each, Eq fl20|) is written as 

dA + (1 - (^)vu - aiu*^v*j = -XiTgA, (21) 
d,u:^-il-a)v*. = -XTgU:^, (22) 
dtUs - (1 - a)vs = -XTgUs- (23) 

The third equation is valid only if strain and rate are colinear, u*j/\us\ = v*j/\vs\- This 
is frequently the case, because any component of Uij not co-linear with Vij relaxes to zero 
quickly. The coefficients a (assuming < a < 1) describes a softening, or more precisely a 
reduced gear ratio: The same shear rate yields a smaller deformation, dtUij = (1 — a)vij ■ ■ ■ , 
but acts also at a smaller stress, aij = (1 — a)7rjj ■ ■ ■ . ai accounts for the fact that shearing 
granular media will change the compression A, implying, as we shall see, dilatancy and 
contractancy. (More Onsager coefficients are permitted by symmetry, but excluded here to 
keep the equations as simple as possible.) 

Next are the continuity equations for density and momentum density, 

dtp + Viipvi) = 0, dtipvi) + Vj{aij + pviVj) = 0, (24) 

where the stress tensor cr.y = P6ij + a*j (with a*j the traceless part) is determined by general 
principles 50. 5lj| as 

P = = (1 - a)PA + Pt - CgVu, (25) 



a, 



^-cy)7i*j-aiu*jPA-%v*j, (26) 



1 - a)7!-s - UiUsPa - Vg^s- (27) 



Again, the third equation (with cr^ = CTijCTij) is valid only if vr^*- and v*j are colinear, 7rj* /|7rs| = 
v*j/\vs\- The pressure P and shear stress as contain elastic contributions ~ tTsjPa from 
Eq (ITUIl . the seismic pressure Pt ~ Tg from Eq (IT^ . and viscous contributions ~ VgXg- 
The off-diagonal Onsager coefficients a, ai (introduced in the equation for the elastic strain 
Uij) soften and mix the elastic stress components. The term preceded by ai is smaller by 
an order in the elastic strain, and may frequently be neglected. 
The balance equation for granular entropy Sg = bpTg is 

dtSg + Vi{SgVi-KViTg)= (28) 

iVgVs + CgVu - 7T^)/Tg- 
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Here, SgVi is the convective, and —nViTg the diffusive ffux. rjgvl accounts for viscous heating, 
for the increase of Tg because macroscopic shear rates jiggle the grains. A compressional 
rate Cc,'^« does the same, though not as efficiently |66|. The term —■jTg accounts for the 
relaxation of Tg, ie., for the conversion of intergranular energy into inner granular one. 



2. Transport Coefficients 

All coefficients a,ai,r]g,(g are functions of the state variables, Uij, Tg and p. As the 
hydrodynamic formalism only delivers the structure of the dynamics, not the functional 
dependence of the transport coefficients, these are to be obtained (same as the energy) from 
experiments, in a trial-and-error process. And the specification below is what we at present 
believe to be the appropriate ones. Generally speaking, we find strain dependence to be 
weak - plausibly so because the strain is a small quantity. One expand in it, keeping only 
the constant terms. We also expand in Tg, but eliminate the constant terms, because we 
assume granular media are fully elastic for Tg — )■ 0, implying the force balance VjCij = pgi 
should reduce to the equilibrium condition, Eq (|3]). Therefore we take a,ai,rig,(g, Kg to 
vanish for Tg — )■ 0. In addition, we also need a, ai to saturate at an elevated Tg, such that 
rate-independence maybe established in the hypoplastic regime. Hence 

Vg = VlTg, Cg = ClTg, K, = KiTg, (29) 

a/a = ai/ai = Tg/{Ta + Tg), 
with a, ai,rji,Qi, Ki, Tq, functions of p only, or the packing fraction 0. Expanding 7 in T^, 

7 = 70 + 71^9, (30) 



astic to the quasi-elastic regime, see Eq ([35 



53| . this does not hold for denser ones. 



we keep 70, because the reason that led to Eqs (1291) does not apply. More importantly, 70 
ensures a smooth transition from the hypop 
below. (Although 70 = in rarefied systems 

The transport coefficients are also functions of p, containing especially a divergent /vanishing 
part ~ (pcp — p)- Assuming that, at p = pcp, the plastic phenomena of stress relaxation, 
softening and dilatancy vanish, Tg relaxes instantly, and the system is infinitely viscous, we 
take 

A, Ai, a, ai, 7f \ r]^^ ~ Pcp - P- (31) 
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We stand behind the temperature dependence with much more confidence than that of the 
density, for two reasons: First, p is not a small quantity that one may expand in, and we lack 
the general arguments employed to extract the T^,— dependence. Second, not coincident ally, 
the p dependence does not appear universal: The above dependence of 71, r/i seems to fit 
glass beads data, whil e 7 1 ~ (pcp — p) "^'^j Vi ~ (Pcp ~ p)^^'^ appear more suitable for 



polystyrene beads, see [55 1. 



At given shear rates, Vs = const, the stationary state of Eq fl28l) - characterized by 
dtSg = 0, with viscous heating balancing T^-relaxation - is quickly arrived at, say within 
10^'^ s in dense granular media, implying 

-fih^r; =vlvi+vlCi, (32) 
where ^ = 1 + -fo/{-fiTg). (33) 

Taking the density for simplicity as either constant or slowly changing, f|f ^ 0, we have a 
quadratic regime for small Tg and low Vg, and a linear one at elevated Tg,Vs- 

Tg = \vsWnihi for 7i7;>7o, (34) 
Tg = v] (r7i/7o) for 7iTg < 70. (35) 

As mentioned above and discussed in the next section, the linear regime is hypoplastic, in 
which the system displays elasto-plastic behavior and the hypoplastic model holds. In the 
quadratic regime, because ~ f ^ ^ is quadratically small, the behavior is quasi-elastic, 
quasi-static, with slow, consecutive visit of static stress distributions. Note that we have 
h = 1 'm the hypoplastic regime, and h 00 m the quasi-static one. 

Eqs (I21|22ll23p also have a stationary solution, 9jA, dtUg = 0, in which the deformation 
rates Vs,vee are compensated by the relaxation ~ Tg. As a result, A = Ac,Us = Uc remain 
constant, and with them also the pressure and shear stress, P = Pc,crs = cfc- This ideally 
plastic behavior is the critical state. In the linear regime, Tg ~ \vs\, both and (Tc are rate- 
independent: A higher shear rate deforms more quickly, but the relaxation ~ Tg increases by 
the same amount. We shall consider the critical state in greater detail in Sec IIV Bl restating 
here only that since the rate-independent critical state is a motion in the linear regime, and 
since it is irreversible and strongly dissipative, it cannot be quasi-static. 
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E. Summary and Three Useful Equations 



With the above set of equations derived, the expressions for energy density and transport 
coefficients in large part specified, GSH is a fairly well-defined theory. It contains clear 
ramifications and provides little leeway for retrospective adaptation to observations. As a 
first step to coming to terms with its ramifications, we examine its basic features. 

Granular rheology as observed may be divided into three shear rate regimes: Bagnold for 
high, hypoplastic for low, and quasi-elastic (ie. quasi-static) for even lower ones. Fast dense 
fiow is in the first regime, in which pressure and shear stress are proportional to shear rate 
squared, p, cr^ ~ f ^. Various elasto-plastic motions, observed especially in triaxial appara- 
tuses, are in the second, rate-independent regime. The third regime is elastic. Static stress 
distribution and elastic waves belong here. This third regime is again rate- independent, 
though the transition between both cannot be. 

Soil mechanics textbooks do not acknowledge the existence of a third rate regime, as they 
take the hypoplastic regime to be quasi-static. As mentioned above and discussed in detail 
in Sec IIII A[ this cannot be right because it is irreversible and strongly dissipative. 

GSH is constructed such that any deviation from elasticity - encapsulated in the coeffi- 
cients a,ai,rig,Q, i^g and the relaxation rate XTg,\iTg of Eqs ( 12T|22|) - vanishes with Tg. 
For Tg = 0, we have dfUij = Vij = dteij, or 

implying perfect elasticity. At very low shear rates, deviations from elasticity are quadrati- 
cally small and negligible, Tg ~ f^. This is the quasi- elastic, or quasi-static regime, because 
the slow motion visits a series of equilibrium, elastic states. 

When Tg is more elevated, we are in the linear regime, Tg ~ \vs\, see Eq f l34p . Here, the 
full complexity of granular media emerges. Nevertheless, true to our starting assumptions 
on granular media: two-stage irreversibility and variable transient elasticity, three scalar 
equations suffice to account for most phenomena. Two account for transient elasticity, 
Eqs fl2T]|23|) . and one for Tg. Eq (|28|). 

dtUs - (1 - (y)vs = -XTgUs, (37) 
dtA + (1 - a)vu - aiUsVs = -XiTgA, (38) 
bpdtTg - K,TgV% = y^.vl - ^,h'{T; - Tl). (39) 
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Some simplifications and one modification were made to arrive at these equations: (1) The 
gradient of Tg is assumed to be small, and linearized in; all other variables were taken to be 
uniform. (2) T^'s convective term is taken to be negligible. (3) An extra source term 7i/i^T^ 
is added to account for an ambient temperature - external perturbations such as given 
by a sound field or by tapping. Generally speaking, any source mechanism contributing 
to Tg is already included in the expression without T^. For instance, given a sound field - 
generated either by loudspeakers or tapping - there is the term on the right hand side of 
Eq ([28]), Cl(^^|^"'^)^ where wf^"""'^ is the fast varying compressional rate of the sound field. 
Coarse- graining it, we may set 

{(Mr')') ^ lih'T', ^ r^^v'^, (40) 

to quantify this contribution, either in terms of the ambient temperature Ta, or a shear rate 
vt needed to produce this T^. Adding this term is a convenient short cut to account for a 
general perturbation, for an ambient temperature without specifying the cause. 

The information on the elastic strain A{t),Us{t) is, for given density, identical as that of 
the elastic stress PA(t),'^sit), because they are always given by the hyper-elastic relation, 
Eq ([I]), or by Eqs ([9l[T0]) for the elastic energy Eq ([8]). The total stress includes the seismic 
pressure Pt and the viscous contributions (of which the compressional one is neglected 
below). We write (assuming that the lowest order terms ~ A^-^ dominate and neglecting 
the term — oiMsPa, of order A^'^) 

P=(l-a)P^ + lT^apH/{p,p-p), (41) 
(T^ = (1 - a)7r, - T]iTgVs. 

[Given Eq ([H]), vr^*-, u*j and a*j are colinear.] In the linear regime, Tg ~ \vs\, the elastic terms 
Pa, are (as discussed at the end of the last section) rate- independent, while Pt ~ T^ and 
rjgVs = TjiTgVs are quadratic in Vs- So both may be written as ei + 62^^, implying a quadratic 
dependence on the rate for 62^^ ^ ei, and rate-independence for e2vl <^ Ci. The first limit 
may be identified with the Bagnold regime, the second with the hypoplastic one. For the 
ease of future references, we note for the hypoplastic regime: 

P = (1 - a)PA, Os = {l- a)'Ks. (42) 

This ends the brief presentation of GSH. 
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III. THE QUASI-ELASTIC REGIME 



A. Quasi-Elastic versus Hypoplastic Regime 

Textbooks on soil mechanics take granular motion in the hypoplastic regime - say the 
approach to the critical state - to be quasi-static. We do not beheve this is right, because 
although slow and rate-independent, it is also strongly dissipative and irreversible. 

Quasi-static motion is never dissipative. Think of sound propagation in any system (such 
as Newtonian liquid, elastic media or liquid crystals), where the velocity is an order in 
the frequency lower than the damping. So sound waves are less damped the smaller the 
frequency is. This is a general feature: Changing a state variable: A slowly, dissipation 
vanishes with the rate of change dfA. The motion is therefore rate-independent in the very 
slow limit, in which dissipation may be neglected. We call it quasi-static because the system 
is at this rate visiting static, equilibrium states consecutively. 

Granular systems are both dissipative and rate-independent in the hypoplastic regime. 
Because of rate-independence, reactive and dissipative terms are of the same order in the 
frequency, and comparable in size. (They are exactly equal in the critical state.) If there 
were only the hypoplastic regime, elastic waves would always be overdamped. Since this is 
not the case, there must be a dissipation-free, quasi-static one that we term quasi-elastic. 

A frequent suggestion is to take a small incremental strain (such as given in an elastic 
wave) to be elastic and free of dissipation, but a large one as elasto-plastic and dissipative. 
We believe this is the wrong way out, because it is (taken literally) illogical and incompatible 
with the notion of a quasi-static motion: Starting from a static state of given stress, and 
applying a small incremental strain that is elastic, the system is again in a static state and 
an equally valid starting point. The next small increment must therefore also be purely 
elastic. Many consecutive small increments yield a large change in strain, and if the small 
ones are not dissipative, neither can their sum be. 

In GSH, it is the strain rate rather than strain amplitude that decides whether the system 
is elastic or elasto-plastic. Of course, small strain increments achieved with a higher but 
short lasting shear rate will indeed provoke elastic responses, if Tg does not have time to get 
to a sufficiently high value to induce plastic responses. Furthermore, the mere existence of 
a quasi-static, quasi-elastic regime does not imply that it is also easily observable. 
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To be specific, we quote a few numbers, tliougli tliese are at best educated guesses. Aside 
from tfie lack of unambiguous experimental data, circumstances are usually complicated by 



density or pressure dependence 46| . We believe, the Bagnold regime starts at shear rates of 
one or two hundred Hz, the hypoplastic regime is say between 10^'^ — IHz, and quasi-elastic 
regime lies possibly below lO^^Hz. 

Finally, we note that backtracing of the stress curve a{t) when reversing the strain, 
e(t) — 7- e(— t), occurs only in the quasi-elastic regime, not the hypoplastic one. (We use a hat 
to indicate a tensor.) The stress is a function of the elastic strain, a = a{u). Reversing u{t) 
will always backtrace d'{t). But only in the quasi-elastic regime may we identify u{t) = e{t). 
Failure to backtrace at hypoplastic rates are not evidence of "history dependence." 

B. A Steep Stress-Strain Trajectory 

As discussed above, in the quadratic regime of very slow shear rates, Tg ~ Iti^p — )■ 0, the 
granular temperature is so small that the system is essential elastic, moving from one elastic, 
equilibrium state to a slightly different one. This is the reason we call it quasi- elastic, or 
quasi-static. Because a — )■ vr and dtu — )■ dti = v, the change of the the shear stress as is well 
approximated by the (hyper-) elastic relation, 

Oum oum duijduke 

Shearing a granular medium at quasi-elastic rates, the result will be a trajectory (T(e) that is 
much steeper than in experiments at hypoplastic rates, such as observed during an approach- 
ing to the critical state. The gradient is given directly by the stiffness constant d'^w/dv?, 
and possibly three to four times as large as the average between loading and unloading at 
hypoplastic rates [because Eq ( 122|) lacks the factor of (1 — a)]. This goes on until the system 
reaches a yield surface of the elastic energy, one of those discussed in Sec III CI We expect the 
system to form shear bands at this point, see Sec IVIII The critical state will not be reached. 
Reversing the shear rate in between will retrace the function o"(t). 

C. Soft Springs versus Step Motors 



Quasi-elastic behavior has not been observed in triaxial apparatus, even at the lowest 
rates. This maybe because they are simply not slow enough. Quite probably though, this is 
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FIG. 3. Why it is hard to observe the quasi-elastic regime if step motors are used, see text. 

also due to the wide usage of step motors in these apphances. Plotting its shear rate versus 
time, 7(t), different shear rates are approximately given as depicted by the two curves of 
Fig 131 Although the curves have different average rates (7), the time-resolved, maximal 
rates 7m are identical. And if the time span of 7m is long enough for Tg to respond, and 7m 
is high enough for the system to be in the linear regime, Tg ~ 7m, the system will display 
consecutive hypoplastic behavior in both cases, irrespective of the average rate (7). 

We suggest two ways here to circumvent this difficulty, both by fixing the stress rate at 
low Tg] As discussed in the last section, a given stress rate has a high shear rate at elevated 
Tg and a low one at vanishing Tg. First is slowly tilting an inclined plane supporting a layer 
of grains. In such a situation, the shear rate remains very small, and the system starts 
flowing only when a yield surface is breached. In contrast, employing a feedback loop in 
a triaxial apparatus to maintain a stress rate would not work well, because the correcting 
motion typically has strain rates that are too high. 

A second method is to insert a very soft spring, even a rubber band, between the granular 
medium and the device moving at a given velocity v to deform it. If the spring is softer by 
a large factor a than the granular medium (which is itself rather soft), it will absorb most 
of the displacement, leaving the granular medium deforming at a rate smaller by the same 
factor a than without the spring. In other words, the soft spring serves as a "stress reservoir" 
for the granular medium. The same physics applies when the feedback loop is connected via 
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a soft spring. Little Tg is then excited, see for instance the experiment discussed in Sec lIVCl 



IV. THE HYPOPLASTIC REGIME 



Hypoplastic motion occurs at an elevated Tg ~ \vs\, in what we have named the linear 
regime. It is rate-independent for given, constant strain rates, in the sense that the increase 
in the stress AcTjj depends only on the increase in the strain, Ae^j = J Vijdt, not how fast 
it takes place. We call this regime hypoplastic because this is where the hypoplastic model 
holds, a state-of-the-art engineering theory that we shall consider in Sec IIVE 1[ 



A. Load and Unload 



In the hypoplastic regime, for given shear rate Vs, the granular temperature relaxes quickly 
to its stationary value Tg = \vs\ ^/vi/li- Inserting this into Eqs (1371 138|) . we arrive at 

dtA = Vs aiUs - \vs\ AiA, (44) 

dtUs = Vs{l - a) - \vs\Aus, (45) 

which are explicitly rate-independent for a = a,ai = di. The last equation is a result 
of inserting the density dependence of Eqs (pTl) and indicates that relaxation of the elastic 
strain becomes slower at higher density, and stops at the close-packed density pcp, where 
the system is essentially elastic. We take A ^ 3.3Ai, as compressional relaxation is typically 



slower than shear relaxation 



66|. 



In this form, it is obvious that loading {vs = \vs\ > 0) and unloading [vs = —\vs\ < 0) 
have different slopes: dtUs/vs = (1 — a) =F i^Us/h). This phenomenon is referred to as 
incremental nonlinearity in soil mechanics, and the reason why no backtracing takes place 
under reversal of shear rate: Starting from isotropic stress, Ug = 0, see Fig HJ the gradient 
is at first (1 — a), becoming smaller as Ug grows, until it is zero, in the stationary case 
dtUg/vg = 0. Unloading now, the slope is (1 — a) + (Aus/h), steeper than it has ever been. 
It is again (1 — a) for Ug = 0, and vanishes for Ug sufficiently negative, see Fig HI Same 
scenario holds for dtA/vg. 

Clearly, only the stress P,(Ts are measurable, not A,^^. The former is calculated em- 
ploying Eq ( 142|) when the latter is given. The resultant expressions can be complicated 
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FIG. 4. The hysteretic change of the shear stress (~ Us) with the strain, as accounted for by 
Eq (j45p . The sign of the shear rate Vs{t) is given in (b), the shear deformation Eg = fQVs{t')dt' in 
(c). Inset (d) is the the temporal evolution of Ug. 

(especially if the pressure is held constant instead of the density, see Sec lIVBj) . but the basic 
physics remains the same - an illustration of why Uij is the better state variable. 

In systematic studies employing discrete numerical simulation, Roux and coworkers have 
obtained great knowledge about the mesoscopic physics on intergranular scales, see eg. jo^]. 
And they were especially able to distinguish between two types of strain, I and II, complete 
with two regimes in which either dominates. However, although type I strain may clearly 
be identified as our state variable Uij, one needs to be aware that regime I is not necessarily 
quasi-static, or quasi-elastic as considered in Sec IIIII The difference is: The relaxation term 
may be temporarily small at hypoplastic shear rates, say because Ug or p^p — p are, see 
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FIG. 5. Loci of static yield surface and the critical states, calculated employing the more general 
energy of Eq (fT8|) . Left: in the space spanned by the three stress eigenvalues, cri, (T2, as; right: in the 
vr-plane of constant pressure, P = CT1+CJ2+CT3, where V^iri = (o"3— cJ2)/P, V^it2 = (2o"i— 1T2— 173)/^. 

Eqs (l44l |46|) . they do not stay small if one wanders in the variable space. At quasi-elastic 
rates, deformation are always free of dissipation. 



B. The Critical State 



1. Stationary Elastic Solution 



When there is complete compensation of the shear rate ~ Vs and the relaxation ~ Tg, 
the stationary solution of Eqs ( l4^E5i) for the elastic strain Uij holds. It is generally called 



the critical state, see [68|, and may be considered ideally plastic, because a shear rate does 
not lead to a stress increase. Setting dtA, dtUg = in Eqs fl38ll37p . we obtain the somewhat 
more general expressions (useful for Sec IIV B 4p , 



Ur. 



aVg — a Ac «i |fs| «! 

Tq A ' \uc\ Ai T„ Ai' 



(47) 



From Eq (122|) . the colinearity of the critical strain and rate, ■u*,|c/|mc| = v*j/\vs\, is easy to 
see. In the hypoplastic regime (for /i = 1, a = a, ai = ai), Uc, depend only on the density 
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and is rate-independent. The critical stress is given by inserting Mc, Ac into Eqs (H2l) . 



PI = Pa ( Ae, Uc) = Va,{BA, + Aul/2A,) , (48) 

TTe = 7Ts ( Ae, Uc) = -2^V^c ^c, (49) 

Pl/vTe = {B/2A)AJuc + uj4A„ (50) 

Pe = (1 - a)Pl, a, = (1 - a)7r,. (51) 



The critical ratio ctc/Pc - same as the Coulomb yield of Eq f|T3|) - is also frequently associated 
with a friction angle. Since one is relevant for vanishing ~ f ^ — )■ 0, while the other requires 
an elevated Tg ~ \vs\, it is appropriate to identify one as the static friction angle, and the 
other as the dynamic one. The dynamic friction angle is always smaller than the static one, 
see FiglHl because the critical state is elastic, and must stay below Coulomb yield, 

Ai/ai < ^/2B/A. (52) 

Textbooks on soil mechanics frequently mention that the friction angle is essentially 
independent of the density - although they do not, as a rule, distinguish between the dynamic 
and the static one, cf. Sec III Ci We assume, for lack of more discriminating information, 
that both are. Therefore, we take ai ~ (pcp — p), because Ai also does, see Eq (H6l) . Quite 
generally, we note that accepting the density dependence of Eqs ( pTl) . we have Ac,Mc being 
monotonically increasing functions of l/(pcp — p)- Same holds for PcCTc ~ -B, though B's 
density dependence make the increase slightly faster. 



2. Approach to the Critical State 

Solving Eqs fl44|45l) for Us(t), A(t), at constant p, Vg and h = l,a = a,ai = di, with the 
initial conditions: A = Ao,Us = 0, the approach to the critical state is given as 

Us{t)=Ucil-e~^'^), es = Vst, (53) 

A(t)= A,(l + /ie-'^^» + /2e-^^^0, 
. _ Ai . _ Ao A 



A-Ai A, A-Ai 

showing that the approach is a simple exponential decay for Us, and a sum of two decays for 
A. It is useful, and demystifying, that a simple, analytical solution in terms of the elastic 
strain exists. Because A ^ 3.3Ai [see the remarks below Eq (H6|) ]. the decay of Ug and /i 
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FIG. 6. Triaxial test curves computed with GSH for a foose (dashed) and dense initial density - 
the axial stress a2 = C3 is being hold constant, with q = ai — (T2- 

are faster than that of f2- Note /2 may be negative, and A(t) is then not monotonic. The 
associated pressure and shear stress are those of Eqs ( 148]I49II50II5T1) . For negative /2, neither 
the pressure nor the shear stress is monotonic. 

For the system to complete the approach to the critical state, the yield surface must not 
be breached during the non-monotonic course of the shear stress. If it happens, shear bands 
will form, see Sec IVIII below, and the uniform critical state will not be reached. 

3. Pressure- Controlled Approach 

Frequently, the critical state is not approached at constant density (ie. volume), but at 
constant pressure P (or a stress eigenvalue Uj). The circumstances are then more compli- 
cated. As A, Ms approach Ac, Uc, the density compensates to maintain P{p, A,Us). Along 
with p, the coefficients a,Q;i,A,Ai, all functions of p, also change with time. In addition, 
with p changing, compressional flow vee = —dtp/p no longer vanishes (though it is still 
small). Analytic solutions do not seem feasible now, but numerical ones are. Reassuringly, 
our result is a perfect textbook illustration, see Fig [61 
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Generally speaking, we have four scalar state variables: p,Tg,Us, A, each with an equation 
of motion that depends on the rates Vg, vu and the variables themselves. In addition, P, as 
are functions of p,Us,A. In the last section, both rates were given, vu = 0, Vg = const. 
As a result, we have p = const, while A(t) and Us{t) were calculated taking the coefficients 
a(p), q;i(p), A(p), Ai(p) as constant. The stress components were then obtained as dependent 
functions. A pressure-controlled experiment means that only the shear rate Vg is given. 
Holding P{p,Us, A) = const (or analogously ui) implies the density p (and with it also 
vee = —dtp/p) is a dependent function, p = p{P,Us, A). Now, in the equations of motion 
for Us and A, one first eliminates vee employing vu = —dtp/p, then eliminates both dtp/p 
and the p-dependence of a(p), ai(p), A(p), Ai(p) employing p = p{P,Us, A). This changes 
the differential equations - which are then solved numerically. 

Many well-known features of Fig |6] can be understood assuming the solutions of Eq ( |53|) 
remain valid, say because the initial density is close to the critical one, hence it does not 
change much during the approach to the critical state. As a result, we may approximate 
a{p), ai{p) , A{p) , Ai(p) as constant, and take vu ~ 0. In addition, we shall assume A ^ Ai, 
or A/(A — Ai) ^ 1, instead of ~ 1.5. Then /2 has the same sign as Aq — A^. Now consider 
the initial stress, cr^ = and P ~ i3(po)AQ^. For P given and B{p) a monotonically 
increasing function of p, the pair Aq — Ac and po — pc have reversed signs. Therefore, we 
have a monotonic change of density for /2 > 0, Aq > A^, po < Pc, and non-monotonic change 
otherwise. At the beginning, the faster relaxation of /i dominates, so A always decreases, 
and p always increases, irrespective of po. After fi has run its course, p goes on increasing 
for the low density case but switches to decreasing for the high-density case, until the critical 
state is reached. These are the cases respectively referred to as contractancy and dilatancy. 
This seems to be indeed what happens, although we do not have A ^ Ai in reality. 

The shear stress cTg ~ o"! — a2 always increases first with m^, until Us is close to Uc- The 
subsequent behavior depends on what A does. With P ~ i3(po)AQ ^ given, as ~ EA^-^ ~ 
P/ A keeps growing if A decreases [loose case, /2 > 0], but becomes smaller again, displaying 
a peak, if A grows [dense case, /2 < 0]. 

As mentioned, during the approach to the critical case, the yield surface may be breached. 
The system will then interrupt its approach to the critical state, and develop shear band 
instead, see Sec I VI II 
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4- External Perturbation or Ambient Temperature 



If one perturbs the system, say by exposing it to a standing sound field, or more simply. 



69| . the critical state is strongly 



by tapping it periodically, such as in a recent experiment 
modified. This is the kind of games/explorations physicists play. To engineers, it may seem 
less than serious, even a bit frivolous, but it does time and again lead to useful insights. 
In 69j, a surprising rate-dependence of the critical shear stress was observed. The stress 
decreases with the tapping amplitude, and increases with the shear rate, such that the 
decrease is compensated at higher rates. Clearly, engineering theories, be they elasto-plastic 
or hypoplastic, that build in rate- independence from the start, could not possibly account 
for this observation. GSH, on the other hand, should be able to, if it indeed provides a 
wide-range qualitative account of granular behavior. 

The consideration of the critical state in the previous three sections takes any granular 
temperature Tg to be a result of the given shear rate, hence Tg = \vs \ ^/vi/lij Eq (IMIl . 
holds. This is no longer the case here, as sound field or tapping will in addition contribute 
to Tg. We have called this portion of Tg the ambient temperature Tq, or the shear rate 
vt = Tq 1^71/771 needed to produce this T^. In Eq (139|) . taking dtTg, ViTg = for a stationary 
and uniform system, we have (for h = 1) 



(54) 



Inserting this into Eq fH7|) . we find the perturbed critical strain and stress: Uc, Ac, (Tc 
■UsV^, given as 




ui Ar ar 1 



ul Ac CTc 1 + v^/ V 



2 u,2- 



(55) 



If there is no tapping, Tq ~ = 0, we retrieve the unperturbed values, Uc = Uc, Ac = 
Ac, (Tc = CTc. With tapping, Uc, Ac, cTc decrease for increasing ~ vt, and increase with 
increasing shear rate Vg, as observed. 
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C. Stress-Controlled Experiments 



1. Long-Lived Temperature and Diverging Strain 

In Sec lIVBt only rate-controlled experiments, with Vg given, were considered. Employing 
Eqs fl38ll37f39l) . we found that 

• the granular temperature quickly becomes a dependent quantity, Tg = \vs\ a/^TiTti, 
with GSH reducing to the hypoplastic model. 

• The exponential relaxation of A,^^ reproduces the approach to the critical state. 

In this section, we examine what happens if we instead hold the shear stress ag constant. 
If Tg = 0, the system stays static, cr^ = const, and there is no dynamics at all. If Tg is 
initially elevated, Ug relaxes, and with it also the stress ag- Maintaining a constant cr^ (or 
similarly, a constant Ug) therefore requires a compensating shear rate Vg. As long as Tg is 
finite, Vg{t) will accumulate, resulting in a growing shear strain es(t) = J Vgdt. As we shall 
see, for Ug close to its critical value Uc, the characteristic time of Tg is ~ (1 — u1/ul)^^ and 
long. Adding in the fact that Tg relaxes algebraically slow rather than exponentially fast, 
the accumulated shear strain can be expected to be very large. 



In a recent experiment, Nguyen et al. [70| pushed the system to a certain shear stress at 
a given and fairly fast rate, producing an elevated Tg. Then, switching to maintaining the 
shear stress, they observed a large total strain eg{t) that appears to diverge logarithmically. 
The authors referred to this phenomenon as creeping, and took it to be a compelling evidence 
that the slow motion of the experiment contains a dynamics and cannot be quasi-static, we 
note that this conclusion sits well with a basic contention of GSH, that what is usually taken 
as quasi-static motion is in fact hypoplastic, with an elevated Tg, see Sec IIII Al 

This experiment may in principle be accounted for by the equations of GSH in the hy- 
poplastic regime, though this - due to its complicated geometry and highly nonuniform stress 
distribution - requires solving a set of nonlinear partial differential equations. Here, we only 
engage in a qualitative consideration of shear-stress controlled experiments in the hypoplas- 
tic regime, assuming uniform variables. Also, we assume for algebraic simplicity that it is 
the elastic shear strain Ug that is being kept constant, not the shear stress ag ~ y/Aug. The 
relevant equations are still Eqs (I38|37|39p . 
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At the beginning, as the strain is being ramped up to Ug employing a constant rate vi, 
the granular temperature acquires the elevated value Tq = {vi/h)\Jrji/'^i. Starting at t = 0, 
Us is being held constant. Setting dtUg = in Eq (|37l) . the shear rate needed to compensate 
the stress relaxation is (see also Eq ( 1471) ) 

- -^^^""^-T.h^M (56) 



L - a Uc]l rji 

(Note that with vt = a/7i / i]i as the shear rate needed to produce Tg at given rates, we 
have the relation Vs/vt = Us/uc-) Inserting Eq ( l56i) into Eqs ( 138)1391) yields 

dt/\ = -AiT,(A - Aoo), Aoo = ^cul/ul (57) 
dtTg = -(7i/M[l - hX- (58) 

These equations may be solved analytically, if the coefficients are constant, which they are 
if the density is. The pressure P{t) ~ A(t) will change with time. This is what we consider 
here. Next, in Sec IIVC"3| we take the pressure as a constant, implying time- dependence of 
density and coefficients. Then, as with the critical state considered in Sec IIVB 3t a general 
solution is possible only by numerical methods. 

The first equation accounts for the relaxation of A, from both below and above A^. The 
relaxation is faster the more elevated Tg is. Writing the second equation as dtTg = —AT^, 
setting h = 1, and employing the initial condition = Tq at t = 0, we obtain the solution 



^0 . _ 7i A u 



,2 



Because of Eq (|56|) , we may rewrite the solution as 

Vo C _ 1 — a 

which implies a slowly growing total strain 

rt 



(60) 



e, - eo = / Vsdt = ln(l + Cvot)/C. (61) 
Jo 

As Vs diminish, it will eventually enter the quasi-elastic regime, ■jih'^Tg — )• 70^^, and the last 
bit of Vs relaxes exponentially. Therefore, e<j does not really diverge. 

Assuming a large rate to ramp up the stress, the initial value for the granular temperature 
To is also large. This will quickly let A be fully relaxed, A = Aoo- Fixing Us = uq is then 
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equal to fixing the shear stress, ctq ~ uq^/A^ = {uq/uc)^/^^- And since the critical shear 
stress is cXc ~ ^/A^Uc, one may rewrite the factor as 

1 - ul/ul = 1 - ao/ac. (62) 

If one chooses to keep (Xs ~ V^Us constant from the beginning, irrespective how far A has 
relaxed, one needs to require dtUs = {us/2A)dtA, resulting in a different proportionality 
Vg ~ Tg to be inserted into the equations of motion. The result should be similar. 



Comparable calculation and analysis were carried out in [70|, using a set of scalar equa- 
tions that may roughly be mapped to the present ones. The quantities Tg, •yi/bp and rji/bp 
(standing for granular temperature, T^-relaxation, and T^-production) were referred to as 
fluidity, aging parameter, and rejuvenation parameter. The present consideration is therefore 
not new, but does provide a treatment embedded in GSH, hence is transparent and unified, 
affording a better founded understanding. 



Since the stress distribution in the experiments of [70| is rather nonuniform, there will 
always be some areas with a shear stress close to Uc. And the system will tend to cave in 
there, resulting in a larger strain accumulation than what the average value for would 



predict. We also note that Cvq, observed to be density-independent in [70|, is correlated 
to the friction angle at high flow rates [see the discussion above Eq (155]) ]. but postpone a 
detailed discussion to a future publication. 

In the experiment, a very soft spring was used to couple the fan and the motor. This 
we believe is an essential reason why this experiment turned out as observed. Usually, 
triaxial apparatus with stiff walls are used. And the correcting rates employed by the 
feedback loop to keep the stress constant are of hypoplastic magnitudes. As a result, much 
Tg is excited, and we have the situation of consecutive constant rates, rather than of that 
of constant stress. The soft spring, as discussed in Sec IIII C] enables quasi-static stress 
correction without exciting much Tg. 



2. Stability above the Critical Shear Stress 

In the last section, we discussed how an initial temperature Tq and its associated shear 
rate vo relax, if the shear stress is hold constant at a value smaller than the critical one, 
as < Uc- After the relaxation has completed, the system is in a static and mechanically 
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stable state. The relaxation is slower the closer as is to dc, becoming infinitely slow for 
as = ac, or equivalently, uq = Uc- Then we have Aqo = Vs = vt, see Eqs (15611571) . with 
especially Tg = Tq remaining a constant, see Eq (13^ . This is indistinguishable from the 
rate-controlled, stationary critical state, which we therefore deduce may be maintained at 
both given rate and stress. In addition, we note that mechanical stability is lost at as = ac 
for an elevated Tg. In contrast, for Tg = 0, a granular assembly maintaining a static shear 
stress exceeding the critical value o"c but not yet breaching the yield condition (as given in 
Sec. Ill Cp is stable, though only precariously so: Any Tg sufficiently large, say caused by 
distant rumbling of the earth, will destabilize it. 

For Uq > Uc, or ag > ac, we have A < in Eq fl59|) . and an initial granular temperature 
Tq will become larger with time, until it diverges. This worsens the instability, and quickens 
the mechanical collapse, but it does not imply the general impossibility for static stress 
values larger than ac, because if the initial Tg is too small, of quasi-static values, the factor 
h = a/1 + '^o/liTg of Eq (133|) is large, and with it also Mc = (1 — a)/A ~ /i, see Eq fH6l) . So 
the factor 1 — Uq/u^ remains positive, and a small Tg will relax rather than explode. 

If the initial Tg is large enough, however, it will indeed grow, and with it also the shear 
rate Vs - though they will not diverge, because the system will leave the rate-independent, 
hypoplastic regime, invalidating the above calculation. The total stress then has a more 
general form including viscous terms, see Eq ( HTl) . and Eq ( l82l) below. Holding it constant 
at a value as > ac will again lead to a stationary state, with constant Tg and Vs- 

Although the shear stress instability for ctq > <^c holds only for stress-controlled experi- 
ments, not rate-controlled ones, the distinction is not always clear-cut in experiments. For 
instance, if a step motor is used, and one has a strain versus time curve such as given by 
Fig [3l than the stress is being hold constant at the plateaus, rendering the stability of the 
uniform system fairly precarious even for a strain-controlled approach. This may well be 
the reason why shear band formation is so frequently observed in the case where the initial 
density is high and the non-monotonic stress trajectory exceeds ac, see FigO 

We emphasize that this aspect of granular behavior comes out of GSH quite naturally, 
without us ever having put it there. It results from the interplay between yield and the crit- 
ical state, or more precisely, between the instability of the elastic energy and the stationary 
solution of the elastic strain. This result does not hinge on any functional dependence of 
energy or transport coefficients, only on the general structure of GSH. 
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It is noteworthy that by including a ambient temperature Ta, as in Sec IIV B 4\ 



d,T, = -{J^/bp)h'[{l - ul/ul)T^ - T„2], (63) 
see Eq ( !39l) . shifts the equihbrium values Tg and Vs relax to, from to 

making its effect more pronounced closer to Uc, or as discussed above, Cc- Note that the 
second expression is algebraically identical to Eq f lS^ . 



3. Stress- and Pressure-Control 

In the last two sections, although the shear stress was controlled, the pressure was not. We 
took vu = 0, a controlled rate, to keep the density constant and the calculation analytical. 
Besides, keeping the shear strain Ug constant is not executable experimentally, and we need to 
find an experimental prescription that is. We therefore revisit the situation - to understand 
what happens if both the pressure and shear stress are kept constant. Though the general 
consideration does not appear analytically viable, one simple solution of a realistic situation 
exists: Rewriting Eqs (I38|37p as 

dtA = ai[usVs - UcVt{A/Ac)] - (1 - a)vu, (65) 
dtUs = (1 - a)[vT{us/uc) - Vs], (66) 

and keeping constant A = Aq, Us = uq, we obtain 

uo Uctti [ul Ao\ 
Vs = vt—, vee = vt- - ^] ■ (6') 

Taking a Uq such that for given Aq the bracket vanishes, and Vu = 0, implies constant density. 
Inserting Vs = vtUq/uc into the balance equation for Tg, Eq fl39|) . we again obtain Eq fl58|) 
with (1621) . The only difference is that there is now a clear prescription for the experiment, 
because constant A,Us,p means that pressure and shear stress are kept constant. So one 
may proceed experimentally by applying an arbitrary pressure, then varying the shear stress 
until the density no longer changes. Tg will then as predicted be long-lived - infinitely so for 
as = <Jc, and exploding for > ac- And no shear stress above (Xc is stable if Tg is elevated. 
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D. Visco-Elastic Behavior 



All visco-elastic systems such as polymers have a characteristic time r that separates 
two frequency ranges: fluid-like behavior for ut <^ 1, and solid-like one for ut ^ 1. Like 
granular media, polymers are transiently elastic, though the transiency is constant and not 
variable. The hydrodynamic theory of polymers 42|-|45l| . defining an elastic strain Uij, and 
employing the equation 

dtU*j - V*j = -U*j/T^e, (68) 

is capable of accounting for many visco-elastic phenomena, including shear-thinning/thicken- 
ing, elongational viscosity, the Cox-Merz rule, and the rod-climbing (or Weifienberg) effect. 

As compared to this equation, the granular version, Eq fl2^ . has an extra coefficient a{Tg), 
and its relaxation time varies as r ~ l/Tg, with Tg a dynamic variable. The first difference is 
not qualitative, as it only accounts for an overall softening. The second difference is crucial, 
because (1) the system is truly elastic when the relaxation time diverges for Tg — )■ 0, and 
may sustain a shear stress statically. (That the static shear stress, as in the case of granular 
media, has a upper limit, is an extra complication not of primary concern here.) (2) The 
relation 1/r ~ Tg ~ gives rise rate-independence, while the constancy of t^^ divides the 
frequency into solid and liquid regimes. 

When there is, in granular media, an ambient temperature much larger than the Tg 
produced by the imposed shear rate Vs, or Ta ^ Tg, vt ^ Vg, then polymers and granular 



media are naturally very similar in their behavior, [vt = a/7i / Tji is the shear rate needed 
to produce the ambient temperature T^, see Eq f l40l ) The ambient temperature may be 
maintained by periodic tapping, or is transported by diffusion from another region of great 
granular activity. The point is, Ta enables relaxation of the elastic strain and stress for a 
static system, and implies a vanishing yield stress. 



1. Creep Motion 

In granular media, one frequently observes a shear band, with a boundary between the 
stationary solid and and the shearing fluid part. More careful experiments (see Komatsu 



et al 



7l|, Crassous et al 



72| . and references therein), however, reveals that the transition 



is not discontinuous, and an exponentially decaying creep motion takes place in the solid. 
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To understand this within the framework of GSH, we first note that granular temperature 
being produced in the fluid region will diffuse into the solid one, and is present there as an 
ambient temperature, enabling stress relaxation. This implies a compensating shear rate if 
the stress is to be maintained. The velocity obtained from integrating the shear rate is the 
observed creep motion. 

The ambient temperature will decay in space, so will the compensating shear rate. Cir- 
cumstances are in fact quite similar to that of Sec II V C| though we need to consider stationary 
but spatially nonuniform states here: The shear rate needed to compensate stress relaxation 
is still as given by Eq f l56|) : instead of Eq fl571) . we take A = Aqo; the balance equation for 
Tg lacks the term dfTg but contains the diffusive current, see Eq (!39l) . 

= (7i/«^i)[l - ul/ul] hX ^ Tg/B\ (69) 

So the decay of Tg ~ exp(— a;/i?) is, for an one-dimensional variation along x, exponential. 
Because of Eq fl56|) . the decay of the shear rate and velocity is also exponential, with the 
same characteristic length 5 ~ 1/ ^Jl — ul/ul = 1/ a/1 — cysl^c- That the decay length 
diverges for Og = should not surprise, because the solid region ceases to exist then. (Note 
that Uci Oc are the critical values at the solid density.) 



2. Nonlocal Fluidization 



An striking phenomenon - discovered by Nichol et al. 73|], and followed up by Reddy 
et al. j?^ in a geometry more amenable to systematic evaluation, is the so-called non-local 
fluidization of granular media. In a vessel containing grains, after a shear band is turned 
on, the medium everywhere, even further away from the band, looses its yield stress, and 
the Archimedes law is observed to hold: A ball that was stuck at whatever height without 
the shear band starts to sink or elevate, depending on its density, until it is equal to the 
granular density. 

The explanation within the framework of GSH is: A ball getting stuck in sand builds 
up an elastic shear stress tTs and the associated elastic strain Us around itself. Without an 
ambient temperature, Ta = 0, this stress holds up the ball's weight if it is not too large, and 
the ball is stationary. But with one, strain and stress relax, and the ball starts to move - a 
striking demonstration of the yield stress vanishing. 
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Once the ball is in motion, the elastic strain rate will quickly become stationary in the 
rest frame of the ball, dtUs = 0, reaching a balance between the deformation rate ~ Vs and 
the relaxation ~ TgUg, or ~ TgUs, see Eq (15^ . Replacing Ug with Vg/Tg in the shear stress 
c"s = '^s{us), we turn the elastic stress surrounding the ball into a viscous one, TTs{vs/Tg). If 
Vg is small enough, we may expand in it, and this viscous stress is linear in the shear flow, 
TTs ~ Vg/Tg. Granular media are then Newtonian fluids, with a viscosity ~ ^/Tg. 

This is the basic GSH-explanation, though a word of caution is appropriate here: Any 
hydrodynamic theory starts from the assumption that the resolution of the theory ("pixel 
size") is small compared to system size, but much larger than any micro- and mesoscopic 
lengths - in the present case, especially the grain diameter d. In [tJ], the diameter of the 
probing rod is only 2d. One may hope that averaging over time and runs will restore the 
macroscopic limit, but this is far from certain when the two scales are that close. 



E. Constitutive Relations 

Granular dynamics is frequently modeled employing the strategy of rational mechanics, 
by postulating a function - of the stress CTjj, strain rate vm-, and density p - such that 
the constitutive relation, dtaij = (tij{aij,Vk£, p) holds (where dt is to be replaced by an 
appropriate objective derivative more generally). It forms, together with the continuity 
equation dtp + Vipvi = 0, momentum conservation, dt{pvi) + Vj(o"jj + pviVj) = 0, a closed 
set of equations for (Tjj, the velocity fj, and the density p (or the void ratio e). Both 
hypoplasticity and barodesy considered below belong to this category. These models yield, 
in circumstances where they hold, a realistic account of the complex elasto-plastic motion, 
providing us with highly condensed and intelligently organized empirical data. This enables 
us to validate GSH and reduce the latitude in specifying the energy and transport coefficients. 

At the same time, one needs to be aware of their drawbacks, especially the more hidden 
ones. First of all is the apparent freedom in fixing €. - constrained only by the data one 
considers, not by energy conservation or entropy production that were crucial in deriving 
GSH. This is what we believe the main reason why there are so many competing engineering 
models. Worse, this liberty explodes when one includes gradient terms. So most models 
refrain from the attempt to account for nonuniform situations, say elastic waves. Second, 
in dispensing with the variables Tg and Wjj, and choosing the shortcut via (t{aij,Vki, p), one 

41 



reduces the model's range of validity and looses the benefit of Ui/s simple behavior: First, 
the models of hypoplasticity and barodesy are valid only for Tg ~ so a Tg that is either 
too small or oscillates too fast will invalidate these models, as will a Tg derived from an 
external source, such as considered in Sec IIVB 4[ Second, as the analytical solution of the 
approach to the critical state in Sec IIV B 21 shows, considering Uij - though it is not directly 
measurable - is a highly simplifying intermediate step. The case for Uij is even stronger, 
when considering proportional paths and the barodesy model. 



1. The Hypoplastic Model 

The hypoplastic model starts from the postulated, rate-independent constitutive relation. 



dtaij = HijkiVke + Kj\Jv1 + evj^, (70) 

n 

where Hij^e, -^ij, e are (fairly involved) functions of the stress and packing fraction [4!]. In- 
cremental nonlinearity as discussed in Sec IIV Al is also part of the postulate. The simulated 
granular response is realistic for deformations at constant rates. 

GSH reduces to the hypoplastic model in the hypoplastic regime, for Tg ~ \vs\, a = 
a,ai = «!, PT,riiTgVij — )■ 0. This is because aij = (1 — a)iTij of Eq (l42l) is then, same 
as Uij, a function of Uij,p, and we may write dt<Jmn = {damn/ duij)dtUij + {damn/ dp)dtp. 
Replacing dtp with the first of Eq (12^ . dtUij with Eq ( l22l) . using Eq (IMl) to eliminate Tg, we 
arrive at an equation with the same structure as Eq (!70|) . Our derived result for Hijk£,Aij 
is different from the postulated engineering expressions, and somewhat simpler, but they 



yield very similar response ellipses, see j66|. (Response ellipses are the strain increments as 
the response of the system, given unit stress increments in all directions starting from an 
arbitrary point in the stress space, or vice versa, stress increments as the response for unit 
strain increments.) 



2. Proportional Paths and Barodesy 

Barodesy is a very recent model, again proposed by Kolymbas [tS^. As compared to 
hypoplasticity, it is more modular and better organized, with different parts in €ij taking 
care of specific aspects of granular deformation, especially that of proportional paths. We 
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take PeP and PaP to denote, respectively, proportional strain and stress path. Their behavior 
is summed up by the Goldscheider rule (gr): (1) A P£:P starting from the stress 0"^ = 
is associated with a PcrP. (2) A PeP starting from aij ^ leads asymptotically to the 
corresponding PcrP obtained when starting at o"jj = 0. (The initial value aij = is a 
mathematical idealization, neither easily realized nor part of the empirical data that went 
into GR. We take it cum grano salis.) 

Explanation: Any constant strain rate Vij is a PeP. In the principal strain axes (ei, £2, ^3), 
a constant Vij means the system moves with a constant rate along its direction, with ei/e2 = 
Vi/v2, £2/^3 = V2/V3 independent of time. What GR states is that there exists an associated 
stress path that is also proportional, also a straight line in the principal stress space, that 
there are pairs of strain and stress path which are linked, and if the initial stress value is 
not on the right line, it will converge onto it. 

If GSH is as claimed a broad-ranged theory on granular behavior, we should be able to 
understand GR with it, which is indeed the case. Given any constant rate Vij, the elastic 
strain will - irrespective of its initial value, relax into the stationary state of Eqs ( I38]l3711 . 

1 - a Ac ai 1 - avu 

Uc = — r— , — = -:- + — 7 , (71) 

A Uc Ai UcAi Vs 

with u*j\c/uc = v*j/vs. Adding in the information from Eqs (|9][T0|) . we also find 

aya,ip)=vyvs. (72) 

If the strain path is isochoric, with vu = and p = const, both the deviatoric strain and 
stress are dots that remain stationary and do not walk down a path as time progresses. 
Clearly, these are simply the ideally plastic, stationary, critical state that we considered in 
Sec lIVBl If vu 7^ with the density p[t] changing accordingly, u*j\c and a*j will walk down 
a straight line along v*j/vs, with a velocity determined, respectively, by Uc{p[t]) and as{p[t]). 

Given an initial strain deviating from that prescribed by Eq ( 17T|) . uq 7^ Uc,Ao 7^ Ac, 
Eqs fl38ll37p clearly state that the deviation will exponentially relax, until they vanish - ie., 
the strain and the associated stress will converge onto the prescribed line. All this is very 
well, but GR states that it is the total stress that walks down a straight line. With 

TTij = Pa (p) + (njPAXj/vs] , (73) 

this fact clearly hinges on (tt^/Pa) ^ a function of A/ug [see Eq (fT^ ] - not depending on the 
density. As long as vu <^ w^, we have A^/uc ~ ai/Ai, a combination that we did assumed is 
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density independent, see Eq (13T!) . partially in anticipation of the fact that the friction angle 
of the critical state, a function of (tTs/P), is independent of the density. And vu/vs must 
indeed remain small to avoid hitting either p^p or pip too quickly. 

In 79|, the results of GSH are compared to that of barodesy, with mostly quantitative 



agreement. (The energy of Eq f|T8|) was employed there. So the results are more realistic.) 
When looking at Cj^, it is easy to grasp that the construction of a constitutive relation is 
only for someone with vast experience about granular media. That we could substitute this 
deep knowledge with the notions of variable transient elasticity, giving rise to a theory just 
as capable of accounting for elasto-plastic motion, is eye-opening. It suggests that sand, in 
its qualitative behavior, may be after all neither overly complicated, nor such a rebel against 
general principles. 



F. Elastic Waves 

That elastic waves propagate in granular media 0, [3] is an important fact, because 
it is an unambiguous proof that granular media possess an elastic regime, and behave as 
elastic media in certain parameter ranges. Experimental exploration of the elastic to plastic 
transition would be equally crucial, and elastic waves remain a useful tool for this purpose. 

There is a wide-spread believe that small, quasi-static increments from any equilibrium 
stress state is elastic, but large ones are plastic. As discussed in Sec lIII A[ this assumption is 
illogical, because a large increment is the sum of small ones. In GSH, the parameter that sets 
the boundary between elastic and plastic regime is the granular temperature Tg. We have 
quasi-elastic regime for vanishing Tg v'^, and the hypoplastic one for elevated Tg Vg- 

A perturbation in the elastic strain or stress propagate as a wave only in the quais-elastic 
regime, while it diffuses in the hypoplastic one. More specifically, we shall derive a telegraph 
equation from GSH, with a quantity ~ Tg taking on the role of the electric resistance |77|. It 
defines a characteristic frequency loq ^ Tg, such that elastic perturbations of the frequency 
u diffuse for w ^ wq, and propagate for w ^ wq. In the quasi-elastic regime, uq ^ 0, and 
all perturbations propagate. In the hypoplastic regime, when Tg gets elevated, so does uq, 
pushing the propagating range to ever higher frequencies. Eventually, the associated wave 
length become comparable to the granular diameter, exceeding gsh's range of validity. 
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To derive the telegraph equation, we start with two basic equations of GSH, Eqs ( !22][2^ . 



pdtVi - (1 - a)V mKimkiuli = 0, (74) 
dtu*^ - (1 - a)v^^ = -\T,u*^, (75) 

where Kimki = —d'^w/duimduki- (For simphcity, we concentrate on shear waves, assuming 
vee = 0.) For Tg — )■ 0, both plastic terms XTgU*j and a ~ are negligibly small, such 
that these two equations reproduce conventional elasticity theory. The variation of wave 
velocities c with stress and density is then easily calculated, because is given by the 



eigenvalues of the matrix KimnjQmQn/ (PQ^) {im is the wave vector). The result |6l| agrees 



well with observations 



60| 



There are two ways to crank up Tg. First is to introduce an ambient temperature, such 
as by tapping or a remote shear band, second is to increase the amplitude of the wave mode, 
because its own shear rate also creates Tg. The granular temperature has a characteristic 
time tt = &p/7i, see Eq (!39|) . that is of order 10~^ s in dense media. For simplicity, we 
assume that the wave mode's frequency is much larger than 1/tt , such that Tg and <y{Tg) 
are essentially constant. This implies 

[d^ + XTgd,) u*^ = '^{1 ~ af X (76) 

Concentrating on one wave mode propagating along x, with Cqs the quasi-elastic velocity and 
u ~ ^tqx-iujt amplitude of the associated eigenvector, we obtain the telegraph equation, 

{d^ + XTgdt) u = il- a)'clVl u ^ c^Vl u. (77) 

(The coefficient a accounts for the fact that granular contacts soften with T^, and the effective 
elastic stiffness decreases by (1 — a)^. In the language of electromagnetism, (1 — a) ~^ is a 
dielectric permeability.) Inserting u ~ e*"^^"""^* into Eq ( !77|) . we find 

c2g2 = ^ iujXTg, (78) 
implying diffusion for the low frequency limit, u ^ AT^, 



72 • 
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and propagation for the high-frequency hmit, u ^ XT^ 



u ~ exp —iu [t^ — ]t^ X 




AT, 



(80) 



(81) 



2c 



The first term in the square bracket accounts for wave propagation, the second a decay 
length 2c/ATg, which is frequency-independent if Tg is an ambient temperature. If Tg is 
produced by the elastic wave itself, it varies between Tg ^ v"^ ^ w^g^ ~ and Tg Vs ^ 
depending on the amplitude, and the decay length is strongly frequency dependent. 

A brief wave pulse, arbitrarily strong, can always propagate through granular media if its 
duration is too brief to excite sufficient Tg for the system to enter the hypoplastic regime. 
The duration must be much smaller than the characteristic time bp/'j of Tg, see Eq ( [39|) . 

V. RAPID DENSE FLOW 

In Chapter [IV]on hypoplastic motion, the seismic pressure Pt and the viscous shear stress 
~ rig were neglected. In this chapter, we consider flows in which they are important, even 
dominant. Including them, we are leaving the rate- independent, hypoplastic regime. Being 
quadratic in the shear rate, the correction come on slowly. This may be the reason rate- 
independence was widely perceived as a basic property of granular media in soil mechanics. 

A. Uniform Dense Flow 

1. Constant Density Experiments 

Starting from Eq f HTj) . we first substitute the unspecific elastic contributions (1 — d)PA, 
(1 — d)7!'s with their critical state expressions, Eqs (!5T]) . because these are the steady state 
values the elastic strain will acquire for constant shear rates, 



P = Pc + 




(82) 



2 Pep- P 



For a stationary temperature, Tg = |t>. 



•s 



I a/^^i /71 [see Eq (IMl) ]. we may abbreviate them as 



P = P, + e,vl 



(83) 
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noting that Pc, c"c, ei, 62, being functions of the density are constant if the density is. Taking 

(84) 



as ffc + e2vl 



P Pc + eivl 

as a friction angle, we have a change between two constant values, from OcjPc for Vs — ?■ 0, to 
62/61 for Vs — 00. Whether the change occurs with f^, or slightly more quickly with f^, is 
here of a fairly subtle difference - though it is of course not in Eq flS^ : That pressure and 
shear stress grow as v1 in the fast limit was already observed by Bagnold [80]. Note that 
since no elastic solution is stable for p < pip, we set Pc, o"c = for any density below p^p, 
where the system is always in the Bagnold regime, P, cTs ~ v"^. 



Both ?7i,7i are believed [5^ to diverge for p pcp. Taking rji ~ (pcp — p)^^'^, 7i ~ 
(Pcp — p)"°'^5 or 61, 62 ~ {pcp — implies first the independence of the high rate friction 
angle, 61/62, from (pcp — p), and second, 

P-P^r^v^J(p^p-p)\ (85) 



There are some experimental evidences for both [55], though the data appear different for 



glass and polystyrene beads 



2. Shear Stress Minimum 



83|. 



If the shear experiment is not executed at given density, but rather at given pressure 
P = Pq, circumstances are more complicated. First, the rate dependence of the friction 
angle is the same as that of the shear stress alone, 

crs/Po = {c^c + e2vl)/Po. (86) 

Second, crucially, the density is a function of the rate: Inverting the first of Eq ( l83ll and 
defining t>^ = (Pq — Pc)/^i = /(p); we have p = f~^{vl). As a result, the rate depen- 
dence in pressure controlled experiments hinges on the density dependence of the transport 
coefficients, which combine to form Pc and 61. 



Third, (Ts and the friction angle observed in [85|, |86|, no longer necessarily 

monotonous functions of the shear rate Vs- In this context, it is important to realize 
that in a nonuniform geometry, keeping the volume constant does not usually maintain 
a constant density. So a non-monotonic relation between cTs and v1 may also happen for 
constant volume, especially if the shear rate is strongly nonuniform. 
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To better understand the last point, consider two uniform subvolumes (instead of the 
continuous non- uniformity of the experiments). They are in contact via a flexible membrane, 
such that their total volume Vi + V2 is a constant. Initially, the total system is uniform, 
with both densities equal, pi = p2, and both shear rates vanishing, 71,72 = 0. Now, if 
72 is cranked up, but 71 remains zero, because Pi(pi,7i) = P2{p2,i2), the density must 
change and the membrane will stretch in one direction, typically with p2 decreasing and pi 
increasing. If system 1 is much larger than 2, the stretching of the membrane will not change 
pi much, as a result, Pi(pi, 71) will remain essentially constant. So will P2 = Pi. As a result, 
the pressure-controlled limit holds. More realistically, if both systems are comparable in size, 
an intermediate case between the pressure- and density-controlled limit will take place. As 
only in the strictly density-controlled limit do we have monotonicity of the shear stress, any 
inhomogeneity in the shear rate may result in non-monotonic behavior of the shear stress. 



3. Comparison to Other Models 



First, we compare GSH to the continuum theory that Boquet et al. [5J] developed to 



account for their experiment, see also [53|. The theory includes the Cauchy stress cTjj, and 
a balance equation for the temperature Tq ~ [see Eq f|TT]) ]. For vu = 0, they are: 

a.,=PS,,-vvl (87) 
with P ~ Tg, rjj'jjK y T( 



Comparing these to the above dense flow expressions of Eqs (13911411) . we find agreement 
except for the fact that the elastic contributions Pc, (Tc are missing. 

Next, we compare GSH to the MiDi constitutive relations. Starting from the postulate 
that granular rheology in dense flows is controlled by the dimensionless parameter of inertial 
number, I ~ 7/v^, Pouliquen et al. distilled two locally applicable constitutive relations 



from experiments and simulations, for the density and the friction angle cTs/P, see 8l|, l82| . 



Pcp-P--I, = • (88) 

Identifying 7 — )■ Vg, the first relation may be combined to form P ~ v1/{pcp — p)^, same as 
Eq ( 185|1 if Pc is neglected. This is to be expected, because the inertial number / (as Savage 
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observed 83|]) does not contain any elastic information. The second expression is similar to 
Eq flH^ . and it does contain elastic contributions, because we may identify /ii = odPc-, and 
= 62/61. The transition between the two friction angles is linear in the reduced shear 
rate /, not quadratic as in GSH, though this is as mentioned a subtle difference, see Fig [71 




V (arbitrary unit) 



FIG. 7. Comparison of the MiDi constitutive relations with GSH. 

More troubling is the claim that Eqs f l88p are valid for pressure controlled experiments, 
not volume controlled ones - implying that the friction angle (Ts/Pq tends to a constant for 
large shear rates. In contrast, GSH contends that it behaves as o-g/Po — ?• e2{p)vl/ Pq, see 
Eq ( 186|) . It is generally a constant only if the rate dependence of 62 (p) cancels that of - 
universally. As mentioned above, density dependence of transport coefficients do vary. So 
we expect the friction angle to display a more diverse behavior, and may even diverge. 



On the other hand, the majority of experiments cited in j46| are inhomogeneous, with 
varying shear rate and density, so controlling the volume does not mean holding the local 
density constant. In princinple, of course, GSH is capable of dealing with these situations 
- though before ons can solve this set of partial differential equations for given boundary 
conditions, we need to have clarified the density dependence of all transport coefficients. 

B. Jamming, or Angle of Repose 

Considering shallow flows on an inclined plane and rotating drums, Aranson and Tsimring 
identified the hysteresis of transition, or the delay between jamming and fluidization, as a 
key feature of granular behavior 8J]. Their theory takes the stress as the sum of two 
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parts, one solid-like, gwij, and a rate-dependent fluid one. q is an order parameter that is 
1 for granular solid, and for dense flow. The authors take the friction angle (f, differently 
than above, as the ratio of the solid stress components, rather than the total ones, and 
postulate a free energy f{g) such that ^ = 1 is unstable for large friction angles, if > ips'-, 
while ^ = is unstable for a small ones, (p < (fR. But both are stable in the intermediate 
region, ips > f > Pr- {fs and ipR are referred to as the angle of stability and repose, 
respectively.) Though this theory does not consider variations in p or Tg, and takes the solid 
stress TTij as an input from some other theory, it provides a pivotal insight: The viability, 
even appropriateness, of using a partially bistable energy to account for the hysteresis. 

Turning now to GSH, we flrst note that for rate-controlled experiments in the hypoplastic 
regime, there is no hysteresis, only uniform and continuous paths to the ideally plastic, 
critical state and back, all given by Eqs f l82|) . with the terms ~ v"^ evident only at higher 
rates. The second type of fluidization takes place either at quasi-static rates or stress- 
controlled. It is discontinuous, and happens because a yield surface is breached - say when 
the ratio tt^/Pa is larger than ^/2A/B, see Eq (fT3ll . This is what we consider now. 

On a plane inclined by a slowly increasing angle p, with y denoting the depth of the 
granular layer on the plane, and x along the slope, we take the stress to be a^x = <^yy = 
azz = P, V^CTxy = CTg, ayz = axz = 0. Integrating VjCiy = gip assuming spacial dependence 
only along y, we flnd axy = gsinp J p{y)dy, ayy = gcosp J p{y)dy, or tanp = (Txy/cTyy = 
as/V^P = tts/v^Pa < a/^Z-B, implying a stability angle ps given by 



Jamming, the reverse transition, is a drop of the shear rate Vs to zero, at given shear 
stress and elevated Tg, which therefore takes place as in Sec IIV Cl And because only stress 
values smaller than ac will jam, and come to a stand still, the angle of repose pn is 



On a plane inclined by a slowly decreasing angle p, the seismic and viscous terms ~ Vg are 
small in the vicinity of pr, and were neglected. Note both ac{p), Pc{p) are functions of p, 
the value of which varying with y is not always flxed. But the ratio o^jPc (a function of 
/S.c/uc = ai/Ai) is density independent, see the discussion below Eq (!52|) . 

The difference between the two angles is clearly a question of Tg, with — )• in the 



tanp^s = \/ AjB. 



(89) 



tanyjR = Ocl^^Pc 



(90) 
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first, and elevated in the second, case. We have ipn < (ps, because the critical state is only 
realized and defined, if the yield surface is not breached in an approach to the critical state. 



VI. COMPACTION 

Compaction - a slow increase of the density at given pressure when the system is per- 
turbed - is a ubiquitous phenomenon in granular media, though not a universal one. For 
instance, the density is observed to both in- and decrease approaching the critical state. 
Within the framework of GSH, compaction stems from the basic fact that the elastic com- 
pression A relaxes in the presence of Tg. Keeping the pressure constant, the density must 
increase to compensate for the diminishing A. This is also the case approaching the critical 
state, though with the difference that, in addition to the relaxation that reduces A, the 
applied shear rate Vs increases it - consider Eq ( l38l) . or dtA — aiUsVs = —\iTgA. In ap- 
proaching the critical state, ~ 7^ 0, A may in- or decrease, depending how large Us, A 
are at any moment, see Sec IIV B"3l 

A. Reversible and Irreversible Compaction 

Consider the pressure P = Pa + Pt assuming vanishing shear rate, Vg = 0, with Pa the 
elastic, and Pt the seismic, contribution, see Eqs fllQ|16ll82p . 




where both B and Pt are monotonically increasing functions of p. So the density must get 
larger when A decreases. It is irreversible because the relaxation of A is. 

At small Tg, the relaxation of A is slow, and the seismic pressure Pt may be neglected. 
This is the limit most soil mechanical experiments are in. Here, only irreversible compaction 
is observed. For larger Tg, the seismic pressure must be included. Then the relaxation of 
A for given Tg increases Pt and decreases Pa , such that P = Pa + Pt remains constant. 
After the relaxation of A has run its course, with Pa — )■ 0, if one modifies Tg but maintains 
P = Pt, the density will change in response, in both directions and hence reversibly. Note 
this is where the fiuid equilibrium condition Eq ([5]), obtained by maximizing the true entropy, 
holds. Consequently, the relaxation of A occurs because it increases the entropy. 




Pa = (1 - «)S(p)Al•^ 
PT = T^{ap'b)/2{p,p-p), 



(91) 
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B. History Dependence versus Hidden Variable 



Changing Tg midway at constant P, with A still finite, will mainly lead to a change in A, 



axation of A, in essence 
and interpreted as a 



because the density responds much more slowly. It disrupts the re 
resetting its initial condition. This phenomenon was observed in 
memory effect. Generally speaking, "memory" is usually a result of hidden variables: When 
the system behaves differently in two cases, although all state variables appear to have the 
same values, we speak of memory-, or history-dependence. But an overlooked variable that 
has different values for the two cases will naturally explain the difference. In the case of 
compaction, the manifest and hidden variables are p and A, respectively. 



C. Tapping and the Edward Entropy 

Numerous experiments have shown that tapping leads to irreversible compaction and 
reversible density change, see the review article 89|. It is usually accounted for by the 
specifically tailored granular statistical mechanics [90] and the Edward entropy SEd, ot some 
variant of it. Substituting the volume V for the energy E, and compactivity X for the 
temperature T, this theory employs dV^ = XdSEd as the basic thermodynamic relation for a 



"mechanically stable agglomerate of infinitely rigid grains at rest" [90|. The entropy Ssd is 
obtained by counting the number of possibilities to package grains stably for a given volume, 
equating it to e'^^''. Compaction is taken as an indication of an increasing S Ed- 
Two reasons prompt us to doubt its appropriateness. First, the number of possibilities to 
arrange grains concerns inter-granular degrees of freedom. These are vastly overwhelmed by 
the much more numerous configurations of the inner-granular degrees of freedom. In other 
words, the Edward entropy SEd is a special case of the granular entropy Sg, and as discussed 



in Sec II D[ we always have Sg <^ S. In equilibrium, where Eq ([^ holds, the entropy S is 
maximal, and macroscopic energy minimal. This is unrelated to the number of possibilities 
to package grains. One would be able to neglect S and concentrate on Sg if these two were 
only weakly connected, if the energy decay from 5*^ to S were exceedingly slow. This is not 
the case. The relaxation of Sg or Tg, via inelastic scattering, is a fast process. 

Second, even assuming a weak coupling between S and Sg, the Edward entropy SEd 
would, as defined, still not be a relevant measure: The actual starting point of the Edward 

52 



entropy is the assumption that Sg does not depend on the energy E, which is always zero 
for infinitely rigid, non-interacting grains at rest, however they are packaged. Taking the 
entropy generally as a function of energy and volume, Sg{E, V), we have, quite generally, 

d^^ = ^d^ + ^dV ^ ^dE + ^dV. 

^ dE dV Tg Tg 

Usually, one keeps the volume constant, d^ = 0, and consider the conventional expression, 
dSg = {l/T)dE. Taking instead E = 0,we have d^^ = {P/T)dV, equivalent to the Edward 
expression dV = {T/P)dSg = XdSg. 

This construction ignores three essential points: First, perturbing the system, allowing it 
to explore the phase space, introduces kinetic energy that one must include. Then clearly, 
E ^ 0. Second, because of the Hertz-like contact between grains, little material is deformed 
at first contact, and the compressibility diverges at vanishing compression. This is a geomet- 
ric fact independent of how rigid the bulk material is. Therefore, infinite rigidity is never 
a realistic limit in granular media, and there is always considerable elastic energy stored 
among grains in mechanically stable agglomerates - even at finite perturbation, as long as 
A is not zero. Finally, Ssd as defined is the granular entropy at vanishing granular motion 
and compression. Its phase space is therefore severely constrained. Generally speaking, each 
classical particle has states in a 6D-space, three for positions and three for the velocities. 
exp(5') is the number of states times the Loschmidt's number; exp{Sg) is the number of 
states times the number of grains, and exp{SEd) is the number of states in 3D space (no 
velocities) times the number of grains. Therefore 

SEd <.Sg<.S. (92) 

Going toward equilibrium, a system searches for the greatest number of states to equally 
redistribute its energy. One bears the burden of proof for the claim that it is sensible for 
the system to neglect S and concentrate on SEd- In contrast, GSH identifies compaction as 
a process taking place at finite Ty and compares the true entropy S of macrostates at that 
Tg. It also accounts for entropy increase, by detailing how macroscopic energy decays into 
granular heat, and how this is converted to true heat. 

Reversible and irreversible compaction as accounted for by GSH is a universal granular 
phenomenon. It occurs at given pressure and Tg, however Tg is created. At the same time, 
numerous experiments show that tapping, though especially efficient, is but one way to 
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achieve compaction, leading to results similar to that of other methods 89|]. So it is natural 
to take the consideration of the last section to hold for tapping as well. This should be true 
for gentle tapping, but the connection to stronger ones warrants further scrutiny. 

Gentle tapping leads to granular jiggling and a small Tg, though one that fluctuates in 
time, with periodic flare-ups. As long as Pt may be neglected, A will relax according to the 
momentary value of Tg, haltingly but monotonically. Since the relaxation is a slow process, 
one could average over many taps to yield a coarse-grained account. Given a granular 
column with a free upper surface in the gravitational field, because a given layer is subject 
to a constant pressure, the density will increase to compensate for the diminishing A. The 
characteristic time of A-relaxation diverges towards the end, and is not a constant, see 91 |. 

Stronger tapping leads to a higher Tg, with A relaxing more quickly. Pj- must now 
be included. Periodically, when all grains are at rest, Pt vanishes, and A is necessarily 
increased to maintain the given pressure. This introduces a non-monotonicity into A(t), 
and raises the question, whether the system, when again at an elevated Tg, will pick up the 
relaxation of A where it was left when the system last crushed to a stop. And why it should 
do so. If it does, we can again take tapping as coarse- grainable, intermittent compaction. 
Then, and only then, does GSH provide an understanding for tapping - though this will be 
a transparent, conventional and demystified one. 



VII. SHEAR BANDS 

A shear band is in its essence the coexistence of static granular solid and uniform dense 
flow. In the first, the grains are deformed and at rest, Tg = 0, with all energy being elastic. 
In the second, the grains jiggle, rattle, move macroscopic distances, with Tg ~ Vg and a 
portion of the energy in Tg. 

The transition from the rate-independent critical state to the Bagnold regime of uniform 
dense flow as the shear rate Vs increases go via two different paths, either gradual and 
uniform, as discussed in Sec IV Aj or discontinuous and nonuniform, via shear bands. 

Approaching the critical state with a high initial density, the evolution of the shear stress 
as is non-monotonic, assuming values temporarily larger than ctc- This is where the system 
has a high probability of breaching an instability, either of the elastic energy at a point on 
the yield surface, as discussed in Sec III C\ or that of Tg, as discussed in Sec IIV C 21 The 
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breaching of the elastic energy will happen with certainty if the system is slowly sheared 
in the quasi-elastic regime. After the breach, the density and elastic strain quickly become 
inhomogeneous, because their fluctuations grow exponentially. This goes on until a stable 
state compatible with the boundary conditions is found again - such as one with a shear 
band, consisting of a low-density fluid region in the shear band, and a high density solid 
region outside. The chaotic transition is difficult to account for, but the stable shear band 
is again simple. 

As we shall see, the shear band has a minimal and constant width at a low shearing 
velocity v. If v is higher, the system's behavior depends on the setup. For given pressure, 
the width £ grows linearly with f, implying a constant rate v/i in the liquid phase. As a 
result, the shear stress, a function of the rate, remains independent of v. This faux rate- 
independence goes on until the band covers the whole system, at which point the quadratic 
rate dependence of uniform dense flow sets in. For given volume, the band width remains 
independent of v, but the shear stress grows quadratically with it. The transition to uniform 
dense flow is again discontinuous. It happens when the shear stress exceeds the critical value 
of the solid density, at which point the solid phase is no longer stable. 

To account for the shear band, we connect the fluid and solid solutions already consid- 
ered employing a set of simple connecting conditions. Denoting the solid and fluid parts 
respectively with the superscripts ^ and ^, the conditions are the equality of the pressure, 
shear stress, and chemical potential. 



[The chemical potential is defined as /i = dw/dp, Eq (j2]). The equality holds because 
otherwise a particle current would fiow across the phase boundary.] All three fields have an 
elastic and a seismic contribution, Eqs ( 125)1261) . With P = Pa + Pt, p = Pa + ^J'T, they are 



tdS r>F S F ,S ,,F 

P = P , a^=cr, , p = p . 



(93) 




in addition to as 



2(1 — a)Aus\/A — ijiTgVg, also with two parts. Denoting the width of 
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the shear band as i, and the velocity difference across the shear band as we assume 

in ffuid: Vs = v/ir^Tg, = A^, uf = u^, (98) 
in sohd: a,Tg,Vs = 0. (99) 

In other words, the elastic strain A and Us have critical values in the F-phase, and ap- 
propriate static values in the S-phase. Strictly speaking, the discontinuities at the S — F 
boundary are in p, A, Ug, but not the shear rate, which decays exponentially in S, as a result 
of Tp-diffusion, see Sec IIVD II We neglect this detail in the qualitative discussion below. 



A. The Fluid Region 

The elastic contribution /iA is a very small quantity: In Pa ~ -BA^'^, a large B compen- 
sates a small A^'^, such that Pa is either much larger than, or comparable to, Pt ~ T^. 
Now, Pt is of the order of Pt/ P, but /xa ~ BA^-^ ~ APa is smaller by the factor A, around 
10^'^ — 10^''. Therefore, as long as Pt ^ APa, we have pt ^ /^a, and p^ = p^ reduces to 
Pt = 0, implying the density in the shear band is fixed, 

p^ = pjil + a). (100) 

Measuring this density therefore yields the value of a [that is important in calibrating the 
energy contribution of granular entropy, see Eq ( TT^ ]. Note that given , the elastic pressure 
Pa is also known, because A-^ = Ac{p^),uf = Uc{p^)- 



1. Given Pressure 

Next, we confine the discussion to the case of given external pressure, 

pex ^pS ^pF 

and given velocity difference v across the shear band. This is an intriguing case, because 
P^ and p^ fix both Tg and the shear rate Tg = v j i = Tg^J'^i/ii]i. Given in addition 
f , the width i of the shear band is also fixed. These are all there is to be known about the 
fluid region. Especially the pressure and the shear stress are given as 

as = a,{p^)-r^^Tgv/L (102) 
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FIG. 8. Shear stress ag as functions of the velocity difference v for given pressure, in a simple-shear 
geometry. The offset gives the elastic contribution, cTcipp)] the parabola is the case without a 
shear band. The thick horizontal line depicts the situation with a shear band, of width which is 
smaller towards left, and equal to the system's width L at the right end. The rate-independence 
of as derived from £ adjusting itself such that Tg ~ v/l remains constant for given pressure. 

Remarkably, the system now displays a faux rate-independence: £ adjusts itself such that 
Tg ^ vjl remains constant for given pressure, independent what v is. The parabola of Fig [8] 
depicts Gg. The offset gives the elastic contributions, Oc- The horizontal line is a result of I 
adjusting. It is indeed easy to mistake a shear band for the uniform, critical state. 

Increasing the velocity v at given pressure alters the width ^, as long as it is smaller 
than the width of the total system L. For larger velocities, the system is again uniform, 
without a solid region. And the consideration of Sec IV Al holds. Until this point, the stress 
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is rate- independent, much longer than without a shear band. 

Given the sohd density (which is fixed by the dynamics, see Sec IVII Bl) and the mass 
per unit length M, mass conservation p^{L — t) + p^i = M determines the total width L 
for given pressure P. 

2. Given Total Volume 

At given total volume L, the behavior is quite different. First, because of mass conserva- 
tion, 

p^{L-i)+p^i = M, (103) 

and because p^ , p^ are given in addition to L, the band width i is fixed, irrespective what 
the velocity v is. As a result, both the shear stress and pressure grow as (v/i)"^ ~ v"^, not 
at all rate- independent. The transition to uniform dense fiow happens when Eq fllOSp is 
violated, for <7c{p^) = = Ccip^) + ViTgVg- 

B. The Solid Region 

The solid region is, in comparison, less fixed. The reason is we have the three connecting 
conditions of Eq ( p3|) . and three quantities to be determined, A"^, wf , p^. Yet, because terms 
are of such different magnitudes in p'^ = p^, it fixes p^ instead of giving a relation between 
p^ and p^ . So it is always satisfied, irrespective what value p^ assumes. Therefore, p^ can 
only be a result of the dynamics: When an instability is breached, the density is changed 
until it gets stuck at some value for p^ , at which the system is again stable. Then of course, 
A^,Ug may be determined for given pressure and shear stress. Nevertheless, we do know 

p^ < p^ and p^ < pc (104) 

must hold. The first inequality can be seen from 

ae(p^) > = a,{p^) + r],TgVs > a,{p^). (105) 

The first greater sign is related to the discussion in Sec IIVC2j the equal sign is one connect- 
ing condition; and the second greater sign is a result of rjiTgVs being positive, in addition 
to the fact that ac is a monotonically increasing function of the density, cf. the discussion 
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below Eq (152|) . The second inequality, < pc, holds for given external pressure P, and 
comes from the following consideration: In the critical state, there is only one free parame- 
ter. Once p is given, Ac, Uc, Pc, o"c also are. Alternatively, one may fix the external pressure 
P = Pc = {1 — a)P^{pc), then Pc{P) is a dependent quantity. In the shear band, because 
the density p^ is already fixed, the elastic pressure Pc{p^) will in general be different from 
the external one, P = Pc{pc), and the difference is taken up by the seismic term Pt- Since 
Pt is always positive, we have Pc{p^) < Pc{pc), implying the external pressure has to be so 
large that p^ < pc holds. Otherwise, shear band cannot exist, and the flow is uniform. 

C. Minimal Band Width 

When the velocity v decreases, the above consideration stops to be valid at some point. 
For instance, p^ is no longer given if Pt ^ APa does not hold. More importantly, the 
width £ will decrease with v (for given pressure) only as long as i is larger than a few grain 
diameter. When v decreases further, £ cannot follow, and will remain at a minimal width. 
To account for this, we need the consideration (as yet quite qualitative) given below. 

The phenomenon of clogging implies that a free surface, if small enough, may be stable 
even when facing downward, with a friction angle of up to 180°, much larger than the coulomb 
yield angle. Similarly, shear bands have a a finite width in the limit of vanishing velocities. 
Both are phenomena not accounted for by GSH as given above. This is connected to the fact 
that hydrodynamic theories are only capable of accounting for spatial variations much larger 
than the correlation length - in the case of GSH especially the grain diameter. Nevertheless, 
there is a tried and proven method of qualitatively accounting for small scale effects such that 
blatant inconsistencies are avoided. (It has been employed eg. for the superfiuid transition 
by including the gradient terms of the order parameter's magnitude {q].) 

In our case, we include higher-order gradient terms ~ (V^Mij)^ in the energy that express 
the extra cost of a nonuniform strain field. A length scale on which elastic strains will 
change is thus introduced. Note non-uniform strain fields necessarily exist at the liquid- 
solid interface, and an infinitely narrow shear band is the result of setting the length of 
strain change to zero. Similarly, a non-uniform strain field of the size of the hole's diameter 
is needed for unclogging. 

Because of momentum conservation, Vjcrjj = 0, the stress stays constant in one- 
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dimensional geometries, even across a shear shear band. Therefore, higher order gradient 
terms in the stress would not do the job. Including higher order gradient terms in strain 
rates would also miss the point that a static inhomogeneity of granular deformation costs 
extra energy and is therefore capped [9^. The lack of a length scale in describing shear 
bands is a known problem in soil mechanics. One popular method to introduce it is by 
adding state variables that account for the couple stress and the Crosserat rotation, see 



eg. 



93| . Including additional variables for the sole purpose of solving our present problem, 
however, does seem unwarranted as it leads to a far more complex theory. 

Starting with an addition to the elastic energy ~ (VfcUjj)^ and introducing the conjugate 
variable (pijk = dw/OVkUij, the elastic stress obtains additional terms of the type 

^k(p^jk ~ V\,,-. (106) 

In a shear band of width i, we therefore expect an additional pressure contribution Pi ~ 1/^^, 
which is to be compared with ~ ~ v'^/i'^. Defining A such that Pt = {Av/i)"^, and 
vq such that Pe = {Avo/iY, where vq is a function of the elastic strain and its difference at 
the interface, the total pressure P = P^ = P^ is 

P = P,{p^) + {A/inv' + v',). (107) 

As long as v is fast enough for i to be larger than, say, 30 grain diameter, f <C f may 
be neglected, and the results of the last sections is recovered. But in the vicinity of a few 
grain diameter, it becomes dominant, and fixes the band width to a value independent of 
V. As shear bands are usually observed to be narrow and rate independent, experiments are 
probably typically in this limit. A constant pressure contribution ~ l/£^ will also stabilize 
a free surface of diameter i that is sufficiently small - a subject that we shall consider 
elsewhere, along with a more quantitative consideration of shear bands. 
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